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Abstract 

Experiments at future colliders will attempt to unveil the origin of electroweak 
symmetry breaking in the TeV range. At these energies the Standard Model 
(SM) predictions have to be known precisely in order to disentangle various 
viable scenarios such as supersymmetry and its manifestations. In particular, 
large logarithmic corrections of the scale ratio -y/i/M, where M denotes the 
gauge boson masses, contribute significantly up to and including the two loop 
level. In this paper we review recent progress in the theoretical understanding of 
the electroweak Sudakov corrections at high energies up to subleading accuracy 
in the SM and the minimal supersymmetric SM (MSSM). We discuss the sym- 
metric part of the SM Lagrangian at high energies yielding the effective theory 
employed in the framework of the infrared evolution equation (IREE) method. 
Applications are presented for important SM and MSSM processes relevant for 
the physics program of future linear colliders including higher order purely elec- 
troweak angular dependent corrections. The size of the higher order subleading 
electroweak corrections is found to change cross sections in the several percent 
regime at TeV energies and their inclusion is thus mandatory for predictions of 
high energy processes at future colliders. 
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1 Introduction 



The Standard Model (SM) of particle physics [1, 2, 3, 4, 5] has enjoyed unprece- 
dented success over the last decades. The discovery of the top quark at the Tevatron 
[6, 7] leaves the Higgs particle [8, 9, 10, 11] as the last undiscovered ingredient to 
complete the SM. While it is possible that the SM remains valid up to energies far 
beyond experimental reaches, most theorists view the SM as an effective theory which 
is embedded in a larger theory usually containing unification of the gauge interactions 
at a high scale Mqut- 

This expectation seems well motivated due to the presence of light neutrino masses 
established at Super-Kamiokande [12] in connection with a seesaw mechanism involving 
Mgut- Also coupling unification in the minimal supersymmetric SM (MSSM) points 
to the existence of a higher scale in nature where the forces unify. 

If, however, the SM is the effective low energy theory of a more complete and unified 
theory at Mqut, the hierarchy problem must be taken seriously. Supersymmetry is 
able to stabilize the quadratic divergences in the scalar sector by canceling these terms 
with the corresponding superpartner loop divergences if, and only if, the superpartner 
mass splittings are not much larger than the weak scale. Another possibility currently 
discussed is that there are large extra dimensions at the TeV scale [13, 14], however, 
such a scenario only trades one problem (the existence of a large scale Mgut) for 
another (the existence of large extra dimensions of the "right" size). 

In any case, while the SM works extremely well, it does not explain electroweak 
symmetry breaking (EWSB). A negative mass squared is introduced by hand in the 
SM, but in the larger theory the reason for EWSB is expected to be dynamical such as 
in typical SUGRA models [15]. While many possible extensions of the SM exist, only 
experiments at future colliders will shed light on the origin of EWSB expected to lie 
in the TeV regime. 

At this point a few general remarks about the usage of the expression EWSB are 
appropriate in order to not be misleading. It has been known for some time now that 
the Higgs mechanism does not lead to a breaking of the local gauge invariance on the 
lattice [16]. In general, all vacuum expectation values (v.e.v.'s) of gauge dependent 
operators (such as (O|0(a;)|O)) can be shown to vanish. As was pointed out in Ref. [17], 
the crucial point about the continuum version in the conventional perturbative formu- 
lation of the Higgs mechanism is not as much the existence of a v.e.v. v, but rather the 
existence of a non-trivial orbit minimizing the Higgs-potential. The apparent breaking 
of the original symmetry by v is due to it being a gauge choice (which always breaks 
the gauge symmetry). In other words, if we were to reformulate the full theory in 
terms of only gauge invariant operators, then no symmetry breaking would be visible 
(but of course new operators would occur describing, for instance, the different masses 
of the electroweak gauge bosons). Since it was also shown in Ref. [17] that the differ- 
ence between the manifestly gauge invariant picture and the conventional perturbative 
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formulation vanishes for observables in the small coupling limit, we prefer to use the 
standard terminology and operators in which the original symmetry is hidden. It is in 
this sense, that we use the expression "broken gauge theory" below. 

The high precision measurements of SLC/LEP have limited the room for exten- 
sions of the SM considerably and in general, they cannot deviate from the SM to a 
large extent without evoking so-called conspiracy effects. It would therefore be very 
desirable to have a Icptonic collider at hand in the future in order to answer questions 
posed by discoveries made at the LHC and possibly the Tevatron. In particular, if 
only a hght Higgs is discovered, say at 115 GeV, then it is mandatory to investigate 
all its properties in detail to experimentally establish the Higgs mechanism including 
a possible reconstruction of the potential and of course of the Yukawa couplings. In 
addition one would have to look for additional heavy Higgs-bosons which could easily 
escape detection at the hadronic machines, but can be discovered at the 77-option at 
TESLA [18, 19, 20] up to masses reaching 80 % of the cm. energy. If any supersym- 
metric particle would be found in addition, it is necessary to clarify and/or test the 
relations between couplings and properties of all new particles in as much detail as 
possible in a complementary way to what would already be known by that time. The 
overall importance of leptonic colliders would thus be to clarify the physics responsible 
for the EWSB which in turn means it must be a high precision machine. 

On the theory side this means that effects at the 1 % level should be under control 
in both the SM as well as all extensions that are viable at that point. The purpose of 
the present work is to summarize the recent activities and results relevant on this level 
of precision from electroweak radiative corrections at energies much larger than the 
gauge boson masses and to apply these corrections to processes relevant to the linear 
collider program. This does not mean that the corrections are negligible for hadronic 
machines, however, for the high precision illustrations we focus here on e'^e" machines 
in the TeV range. 

At the expected level of precision required to disentangle new physics effects from 
the SM in the O {< 1%) regime, higher order electroweak radiative corrections cannot 
be ignored at energies in the TeV range. As a consequence, there has been a lot of 
interest recently in the high energy behavior of the SM [21, 22, 23, 24, 25, 26, 27, 28, 29]. 
The largest contribution is contained in electroweak double logarithms (DL) of the 
Sudakov type and a comprehensive treatment of those corrections is given in Ref. [28] 
to all orders. The effects of the mass-gap between the photon and Z-boson has been 
considered in recent publications [30, 31] since spontaneously broken gauge theories 
lead to the exchange of massive gauge bosons. In general one expects the SM to be in 
the unbroken phase at high energies. There are, however, some important differences of 
the electroweak theory with respect to an unbroken gauge theory. Since the physical 
cutoff of the massive gauge bosons is the weak scale M = Mw ~ Mz ~ Mh, pure 
virtual corrections lead to physical cross sections depending on the infrared "cutoff" . 
Only the photon needs to be treated in a semi-inclusive way. Additional complications 
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arise due to the mixing involved to make the mass eigenstates and the fact that at high 
energies, the longitudinal degrees of freedom are not suppressed. Furthermore, since 
the asymptotic states are not group singlets, it is expected that fully inclusive cross 
sections contain Bloch-Nordsieck violating electroweak corrections [33]. 

It has by now been established that the exponentiation of the electroweak Sudakov 
DL calculated in Ref. [28] via the infrared evolution equation method (IREE) [34, 
35] with the fields of the unbroken phase is indeed reproduced by explicit two loop 
calculations with the physical SM fields [30, 31, 36]. One also understands now the 
origin of previous disagreements. The results of Ref. [26], based on fully inclusive 
cross sections in the photon, is not gauge invariant as already pointed out in Ref. 
[28]. The factorization used in Ref. [25] is based on QCD and effectively only takes 
into account contributions from ladder diagrams. In the electroweak theory, the three 
boson vertices, however, do not simply cancel the corresponding group factors of the 
crossed ladder diagrams (as is the case in QCD) and thus, infrared singular terms 
survive for left handed fermions (right handed ones are effectively Abelian) in the 
calculation of Ref. [25]. The IREE method does not encounter any such problems 
since all contributing diagrams are automatically taken into account by determining 
the kernel of the equation in the effective regime above and below the weak scale M. 
It is then possible to calculate corrections in the effective high energy theory in each 
case yielding the same result as calculations in the physical basis. Thus, the mass gap 
between the Z-boson and the photon can be included in a natural way with proper 
matching conditions at the scale M. For longitudinally polarized gauge bosons it was 
shown in Ref. [37, 38] that the leading and subleading (SL) kernel can be obtained 
from the Goldstone boson equivalence theorem. 

We specify in the next section how the high energy effective theory is obtained 
from the SM and illustrate the approach followed in the main part of this work. 

1.1 The Standard Model 

The complete classical Lagrangian £ciass of the electroweak SM (EWSM) reads in 
terms of the physical fields, i.e. the mass and charge eigenstates A^, Z^, VF^, H, /, z/, 
M, and (i, the would-be Goldstone fields 0^ and %, and the physical parameters e, Mw, 
Mz, Mh, m/, and V, as follows [39]: 

f=l,u,u,d i 
f=l,u,u,d i *wCw 
i,j V ^Sw 
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The quantization of the EWSM requires the introduction of a gauge-fixing term 
and of Faddeev-Popov fields. We introduce a gauge-fixing term of the form 



A\2 
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(2) 



with hnear gauge-fixing operators 



(3) 



This general hnear gauge contains five independent gauge parameters ^a, a — A,Z, ±, 



and a = Z, ±, where = ^ 
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For — and = the terms involving the would-be Goldstone fields in (3) 
cancel the mixing terms V^9^0 in the classical Lagrangian (1) up to irrelevant total 
derivatives. This gauge is called 't Hooft-Feynman gauge and is used in the following 
if not stated otherwise. 

The corresponding Faddeev-Popov ghost-field Lagrangian reads 



ie{d''u-^)W^^ 



u 



--w Z 

— u 



- {u+ 



1 



2s, 



iH + ix)u+- 



ct - s 



u — 



W-, 0+ ^ r, i 



(4) 

Adding up all terms (1), (2) and (4) we obtain the complete Lagrangian of the 
EWSM suitable for higher-order calculations, 



'GSW 



ghost 



(5) 



The Feynman rules which can be derived from the Lagrangian in Eq. (5) are given in 
appendix 6.2 in the 't Hooft-Feynman gauge for the physical fields of the broken gauge 
theory. 

At high energies and for processes that are not mass suppressed or dominated 
by resonances we can neglect particle masses and terms connected to the vacuum 
expectation value (v.e.v.) of the broken gauge theory to the level of SL accuracy [40]. 
Thus, instead of the Lagrangian in Eq. (5) we use a high energy approximation of 
^symm which is bascd on the fields of the unbroken phase in the symmetric basis and 
neglect all terms with a mass dimension. It is composed of a Yang- Mills part, a Higgs 
and a fermion part which are given by [41]: 

Cym = -\ {d,W: - d^W; + ge^^^^W'^W^)^ - \ {d,B, - d,B,f (6) 

where e"^'^ is the totally antisymmetric tensor of SU(2). The Higgs part consists of a 
single complex scalar SU(2) doublet field with hypercharge Y — 1: 



with (j)^{x) — {H{x) + and where the v.e.v. is neglected. 0"*", and x 

denote the would-be Goldstone bosons and H the physical Higgs field. They couple to 
the gauge fields via 



£h = iD^<^y {D>^<l>) (8) 

where we omit the self coupling part (the potential) and the covariant derivative is 
given by 

D, = d,-tgT'^W^ + tg'^B^ (9) 

The left handed fermions transform as doublets and the right handed ones as singlets 
under the gauge group. The fermionic part of the symmetric Lagrangian is then given 
by 

>Cf = E {Lj^I'D.L,) + E {R,ii'D,R,) - E {Lfi.iRi^ + h.c) (10) 

3 j j,l 

The covariant derivative acting on right handed fields contains no term proportional 
to g. The Yukawa coupling matrices are denoted by Gji noting that for up-quarks, the 
charge conjugated Higgs field must be used. The high energy effective symmetric part 
of the Lagrangian is then given by 

'^symm = -^YM + + + C^x + -^FP (H) 

where the corresponding ghost and gauge fixing terms are given by 

1 



.-fix 



(12) 



with 
and 



iw 



--dmi , = -^d^'B^ (13) 



>Cfp = TI"(x)^k^(x) (14) 
where is the variation of the gauge fixing operators under the infinitesimal 



gauge transformations characterized by 6'^{x). The Faddeev Popov ghosts are denoted 
by u'^{x). The corresponding Feynman rules are thus analogous to a theory with an 
unbroken SUl{2) x C/y(l) and fermions or scalars in the fundamental representation 
respectively. The new ingredient in Cp is the Yukawa term. In addition we have in the 
gauge boson sector the coupling of the gauge bosons to scalars through the covariant 
derivative in £h- This effective regime corresponds to region 1) in Fig. 1 where the 
wavy line separates the transverse sector (analogous to an unbroken gauge theory) and 
the scalar sector, where for the would-be Goldstone bosons the equivalence theorem 
(E.T.) must be used. Note that all gauge bosons contained in £symm are massless and 
an infrared cutoff will treat PV^ and fields in the same way. 
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In the following we always use the 't Hooft-Feynman gauge. For the high energy 
regime, the Lagrangian in Eq. (11) is convenient since it allows for an approach via 
the IREE method [34, 35] described in section 3 and thus, for a consistent treatment of 
higher order SL corrections at high energies. In Ref. [40] it was proven that Eq. (11) 
at one loop to SL accuracy gives the same results as calculations based on the physical 
Lagrangian in Eq. (5). The approach in Ref. [40] uses coUinear Ward identities to show 
that SL contributions from the v.e.v. part of the Lagrangian (5) do not contribute 
additional terms not already contained in £symm- In particular this means that for 
longitudinal degrees of freedom at high energy we can employ the Goldstone boson 
equivalence theorem and to SL accuracy, we treat the would-be Goldstone bosons 0"*", 
0~ and X as physical degrees of freedom in the ultrarelativistic limit. At higher orders, 
all terms related to the renormalization of the Goldstone bosons are sub-subleading 
(SSL). 

Fig. 1 also indicates that this approach is only valid in the high energy regime 
^/s M and that the QED corrections from below the weak scale must be included 
by appropriate matching conditions at M. 

Thus the overall approach consists of identifying the relevant degrees of freedom 
in region I) and II), integrating out the contributions to SL accuracy and by matching 
the solution found in II in such a way that at the weak scale M the solution in region 
I) is reproduced. 

1.2 Organization of the paper 

The paper is organized as follows. In section 2 we summarize the various ingredients 

needed to calculate SL virtual corrections in unbroken gauge theories. While these 
corrections do not lead to physical observables in those theories, the IREE approach 
allows for an application of the results of section 2 to broken gauge theories in the 
high energy limit in section 3. The QED effects from the region below the weak scale 
are implemented with the appropriate matching conditions as indicated above. As 
mentioned above, we use the term "broken gauge theories" in the sense that the local 
symmetry is hidden due to the degeneracy of the vacuum ground state and thus not 
evident in the physical states. The associated local BRST relations, however, still hold 
[16, 17]. 

In section 4 the results summarized in section 3 are applied to specific processes 
relevant to a future linear collider program. In particular the importance of the higher 
than one loop corrections is emphasized. We present our concluding remarks in section 
5 and discuss lines of future work needed for precision prediction at future TeV colliders. 
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Figure 1: The schematic depiction of the effective high energy regimes {^/s » M » A) 
in the framework of the infrared evolution equation method. In region I), the high 
energy corrections are obtained effectively in the unbroken SUl{2) X Uyi^) theory 
described by £symm in E<1- (H) where all terms connected to the v.e.v. can be neglected 
to SL accuracy. For external fermions and transverse gauge bosons this picture contains 
at the subleading level Yukawa interactions and contributions from Higgs doublets 
to the anomalous scaling violations. For external longitudinal gauge bosons {(f) — 
{4>~^ , 4>~ , x}) the equivalence theorem (E.T.) is employed yielding effectively a scalar 
theory charged under the unbroken gauge group. Again this scenario contains at the 
subleading level Yukawa terms introduced by the spontaneous symmetry breaking. 
For all charged particles, the soft photon effects, regulated here by a fictitious photon 
mass A, are included by integrating in region II) which incorporates pure QED effects 
including mass terms. In the calculation A is replaced by a cutoff /i on the exchanged 
The matching condition is given by the requirement that the high energy solution 
in region I) is obtained if the infrared cutoff is chosen to be the gauge boson mass 
M. 
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2 Unbroken gauge theories 



In this section we summarize the results obtained for virtual corrections in unbroken 
gauge theories at high energies. These contributions will be crucial for the high energy 
regime of the SM in section 3. 

2.1 Sudakov double logarithms 

The high energy asymptotics of electromagnetic processes was calculated many 
years ago within the framework of QED [42]. In particular the amplitude for e+e^ 
elastic scattering at a fixed angle (s ~ |t| ~ |m| S> 3> A^, where m is the electron 
and A a fictitious ^ photon mass) in the DL approximation has the form 

M = M^r^[^,^) (15) 

where A^Born IS the Born amplitude for e+e" scattering and F is the Sudakov form 
factor. The DL approximation applies in the energy regime 

9 9 9 

e , 2 * e , s , m , , 

log ~ 7Z?log-^log— ~ 1 (16) 



47r2 ^ m2 47r2 * A^ 

where the QED coupling e^/An <S 1. Thus each charged external particle effectively 
contributes vT to the total amplitude. The Sudakov form factor appears in the elastic 
scattering of an electron off an external field [42]. It is of the form: 

To specify R it is convenient to use the Sudakov parametrization of the momentum of 
the exchanged virtual photon : 



k = vpi + up2 + k± (18) 



for massless fermions and 



CTJl^ \ I TP? \ 
Pi —P2j+uip2 ^Pij+A;_L (19) 

for massive fermions. pi and p2 are the initial and final momenta of the scattered 
electron and in the following we denote the Euclidean component 



-'^A plays the role of the infrared cut-off. In physical cross sections the divergence in A of the elastic 
amplitude is canceled with the analogous divergences in processes with soft photon emissions. 
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Figure 2: The virtual Sudakov DL-phase space in massive QED for the function R in 
the {u,v} representation. The cutoff /i plays the role of A for /i -C m. The shaded 
area is the region of integration and is symmetric with respect to u and v. For fi > m 
the relevant phase space is mass independent. 

R{^-:^i ^) can then be written as the integral over u and v after rewriting the measure 
as df'k = (Pk±(Pk\\ with 

d^k^ = \k^\d\k±\d(f) = ]^dk\d(p = Tidk\ (21) 
d^k\^ =\d{k\k'')/d{u,v)\dudv^^-dudv (22) 

where we turn the coordinate system such that the pi,p2 plane corresponds to 0, x and 
the y,z coordinates to the k± direction so that it is purely spacelike (see Eq. (20)). 
The last equation follows from pf = 0, i.e. pf^ ~ pf^ and 

Integrating according to the DL phase space of Fig. 2 (where /i plays the role of A): 

r{—,^%] = f'du fdv f ^ \ , ) i ^ \ , ] e{suv-\') 
A'' J Jo Jo \u + m'^v/sj \v + m''u/sj 



12 



^ C — e(u-m^v/s) e(v-m^uls) e{suv - X^) (24) 

Jo u Jo V ^ ^ ^ ' 

where s ~ \t\ ~ 1'P\V2- The first two factors in the integrand correspond to the propa- 
gators of the virtual fermions which occur in the one-loop triangle Sudakov diagram. 
The Q - function appears as a result of the integration of the propagator of the photon 
over its transverse momentum A;^: 

writing it in form of the real and imaginary parts (the principle value is indicated 
by V\ The latter does not contribute to the DL asymptotics and at higher orders 
gives subsubleading contributions. We note that the main contribution comes from 
the region near the photon mass shell: 

SUV = + fci . (26) 

To DL accuracy Eq. (24) gives for A <^ m: 

^ / s rf?\ 1 , a s , s , rr? 
A^ y 2 A^ 

where the result comes equally from two different kinematical regions, v ^ u and 
M ^ V as is evident from Fig. 2. Therefore one can write R — 2r. 

We can obtain physical insight by presenting the two equal contributions separately. 
In the first region, with v ^ u, the virtual photon is emitted along pi and the parameter 
V is given by the ratio of energies of the photon and the initial electron. Here instead 
of u, it is convenient to use Eq. (26) to replace it by the square of the transverse 
momentum component of the photon. Then integrating over v and fc^ according to 
the DL phase space in Fig. 3 gives 

/•! dv r*^^ dk]_ dk\ r"^H\k±\/m,i) 

A/V5 V A2 k]_ + w?v'^ ~ A2 k\ J\kj_\/v^ V 

in the DL approximation, which may be evaluated to give half of R. The quantity r is 
proportional to the probability Wi of the emission of a soft and almost coUinear photon 
from an external particle with energy ^/s and mass rrii, i.e. 

If several charged particles participate in a process, for example e~^e~ — > ffff, then 
analogous contributions appear for each external line, provided all external invariants 
are large and of the same order. This leads to the general result 

M = Msorn exp 1^-^ f2 ^')) (30) 
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Virtual Sudakov DL-Phase Space 




s 



Figure 3: The virtual Sudakov DL-phase space in massive QED for the function R in 
the representation. The shaded area is the region of integration. For fx <^ m 

the cufoff plays the role of A in the text. For fi > m the relevant phase space is mass 
independent as in the {u,v} representation of Fig. 2. 

where n is the number of external lines corresponding to charged particles. In sum- 
mary the soft emissions described by the Sudakov form factor is a quasi-classical effect 
which does not depend on the hard dynamics of the process. In particular there are 
no quantum mechanical interference effects in the DL Sudakov corrections, for large 
scattering angles. 

2.2 Gribov's factorization theorem 

In this section we discuss a factorization theorem due to Gribov [43, 44, 45, 46]. 
It was originally derived in bremsstrahlung off hadrons at high energies in the context 
of QED but can appropriately be extended to non-Abelian gauge theories. We follow 
the original derivation for real emission processes noting that the form of factorization 
of virtual corrections must be analogous due to the KLN theorem [47, 48]. 

Consider bremsstrahlung off a fermion with mass m in the laboratory system. We 
denote the invariants according to the notation depicted in Fig. 4 as follows: 
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Figure 4: Bremsstrahlung in a process involving charged fermions. At high energies 
only the external legs contribute to DL accuracy. 



s'^ip[+p',r, h^{p',-p^Y (31) 

The usual eikonal argument is that for 2pik <^ rn?, only the diagram on the l.h.s. is 
large, yielding (neglecting m in M): 



M^-^M{s',t2) (32) 
Pik 



At large energies we have 



Ifcl 

2piA; = 2|pi| |fc|(l -cos^) + (33) 

|Pi| 

Thus for Ai to be large we need in any case: <^ 1. It follows that the condition 
|pi| 1^16*^ <C should be fulfilled for small emission angles 9. 

Gribov observed, however, that A4 is large in broader region! In the region 2pik > 

m^, < 1 we have with pi = (p?,0,0, |pi|), = (/c°, A;° sin ^, 0, cos^) and e* = 
(0, cos^, ±i, — sin^): 

M^-j^M{s',t2) (34) 

i.e. a sufficient condition is: \k±\ ^ \k\9 <^ m. 

The proof proceeds as follows: 
M is taken on the mass shell in order to ensure gauge invariance! In covariant form 
the conditions read^: 

2P1/C ^ , '2p2k 2 4.pikp2k 2 /grX 

< 1, < 1, K_L = < m (35) 



^We consider here only the case of initial state radiation in analogy to Ref. [43, 44]. 
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We can then write the amphtude in a gauge invariant form: 

£ P2 ^ V\ 



M 



M{s,t) 



\P2k pik 
and write it in the following way: 

Gauge invariance yields: 

PikMi +p2kM2 + qkMs = 



(36) 



(37) 



(38) 



At high energies t2 - h ^ 2qk = ^2qp2 + ^2qpi + 2|fe_L||q^| ^ 2|fe_L||qf^|. Thus, 
2qk is small, so that there is no need to distinguish between s, s' and between ti, t2 in 
the inner on-shell amphtude M of Eq. (36). Eq. (38) then reads 



Thus 



Ml = 
M2 = 



PikMi = -p2kM2 

M{s,t) 



Pik 
M{s,t) 
P2k 



-P2k T{s,t,pik,p2k) 
-pik T{s,t,pik,p2k) 



(39) 

(40) 
(41) 



where the only pole of Mi is at pik, that of M2 at p2k. M3 and T have no singularity 
at pik or p2k, and 



P2c 



-Mm = ( ^ - [M{s,t) +pikp2kT] + q^Ms 



\k 



M{s,t) + -^klsT 



+ e*qM3 



(42) 
(43) 



Since M3 and T are functions of pik, p2k they could be of order of the first term. 

To show that this is not the case, consider the imaginary part (discontinuity) of 
M. in pik, p2k: 



ImM^ --\ki_\s ImT + e*q Im M3 



(44) 



It is determined by all possible splittings (not in s, ti, ^2) like the ones depicted in Fig. 
5. The simplest two particle intermediate state contains the amplitude 



e Ki e pi 
Kik pik 



M,((Pl+^)^/) 



(45) 
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Figure 5: Higher order splittings determining the imaginary part of the scattering 
amphtude. At high energies, Gribov proofed that the pole terms in the variables 2pik 
do not dominate the amplitude and that the large terms factorize with respect to 



where Ki is the momentum of the charged particle. 

As pi is large and (pi + <^ s, |ki| is also large and along p^. 

Pi , Q± 



Kik pik Kik 
where is the component of Ki perpendicular to p^. Thus 



(46) 



M, = ^M,((pl + A;)^g'') (47) 

We therefore observe that the large terms cancel! For higher splittings the cancella- 
tion proceeds analogously since in the intermediate states all particles formed at high 
energies are parallel to the original particle momentum. 

Thus, the large contributions to the original bremsstrahlung amplitude are given 

by 

M = -j^M{s,t) (48) 

for 2pik <ti s, 2p2k <S s and kx^ <tC A* < rn^. The cutoff fi is introduced for later 
convenience. An analogous factorization in then holds for virtual corrections with 

p^ since the sum of real and virtual corrections must be independent of the 

infrared cutoff. 

In order to treat Non-Abelian gauge theories we need to introduce a gauge invariant 
cutoff on all virtual particles with momentum Ki. 

i^i±^ > A*^, i^i = ViPi + UiP2 + Ki_i_ (49) 
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and it is understood that // ^ Aqcd in order to remain in the perturbative regime. 
The crucial point now is that fi determines both the positions of the thresholds in the 
variables 2pik, 2p2k and the minimum momentum transfers [49]: 

+ -) + '^2±\1 + -)- > V (50) 

for <^ /i^. Thus the cut starts from 4//^ and 

\{k - qif\ = {qis. - k±f - s{ui - u){vi - v) > (51) 

Now the same dispersive arguments are applicable to QCD as they were in QED. 

Thus we can consider again the simplest situation, when the additional soft gauge 
boson is emitted in the process with all invariants Sij large. Of course, for the emission 
of a boson almost coUinear to the particle the direction of the particle with momentum 
Pi, the invariant 2kpi is small in comparison with s. In the case of non-Abelian gauge 
theories the corresponding amplitude for the emission of a soft gauge boson with small 
k'^j_ <^ ji^ has, according to the Gribov theorem as derived above, the following form 
in non-Abelian theories: 

M\p,,...,Pn,k■^Ji') = j29s^^T\J)M{p,,...,Pn.^Ji^). (52) 

where gs denotes the QCD (or SU{N)) coupling. The possible corrections to this 
factorized expression are of the order of k\/iJi'^. However, to DL accuracy, we can 
substitute /x^ in the arguments of the scattering amplitudes by its boundary value 
k\. Notice that the amplitude on the r.h.s. of (52) is taken on-the-mass shell, which 
guarantees its gauge invariance. The result (52) is highly non-trivial in the Feynman 
diagram approach. It means, that the region of applicability of the classical formulas 
for the Bremsstrahlung amplitudes is significantly enlarged at high energies. 

The form of the virtual factorization and the subsequent resummation is the topic 
of the following section. 

2.3 Infrared evolution equations 

Sudakov effects have been widely discussed for non-Abelian gauge theories, such 
as SU{N) and can be calculated in various ways (see, for instance, [50, 51, 52, 53, 54, 
55, 56, 57, 58]). We consider here the scattering amplitude in the simplest kinematics 
when all its invariants Sij = 2piPj are large and of the same order Sij ~ s. A general 
method of finding the DL asymptotics (not only of the Sudakov type) is based on 
the infrared evolution equations describing the dependence of the amplitudes on the 
infrared cutoff /i of the virtual particle transverse momenta [34, 35] . This cutoff plays 
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the same role as A in QED, but, unlike A, it is not necessary that it vanishes and 
it may take an arbitrary value. It can be introduced in a gauge invariant way by 
working, for instance, in a finite phase space volume in the transverse direction with 
linear size / ~ l/ii. Instead of calculating asymptotics of particular Feynman diagrams 
and summing these asymptotics for a process with n external lines it is convenient to 
extract the virtual particle with the smallest value of \kj_\ {k± ± Pj,Pi) in such a way, 
that the transverse momenta \k'j_\ of the other virtual particles are much bigger 

k'l > fei > //^ (53) 

For the other particles k']_ plays the role of the initial infrared cut-off jj,^. 

In particular, the Sudakov DL corrections are related to the exchange of soft gauge 
bosons, see Fig. 1. For this case the integral over the momentum k of the soft (i.e. 
|/c°| -C ^/s) virtual boson with the smallest k± can be factored off, which leads to the 
following infrared evolution equation: 

2 (27r)4 Js^kl^^^ + ie {kpj){kpi) 

xr"(j)nOM(pi,...,P„;fel), (54) 

where the amplitude A4(pi, ^1) on the right hand side is to be taken on the mass 

shell, but with the substituted infrared cutoff: /i^ — > k\. From Eq. (18) and the 
on-shell condition (26) it is clear that (^^pi'^'^.p^^ — and that Eq. (54) has the required 

factorized form for the virtual corrections according to the discussion in section 2.2. 

The generator T"(Z)(a = 1, ...,A^) acts on the color indices of the particle with 
momentum pi. The non-Abelian gauge coupling is g. In Eq. (54), and below, k± 
denotes the component of the gauge boson momentum k transverse to the particle 
emitting this boson. Note that in Sudakov DL corrections there are no interference 
effects, so that we can talk about the emission (and absorption) of a gauge boson by 
a definite (external) particle, namely by a particle with momentum almost collinear 
to k. It can be expressed in invariant form as fe^ = mm{2{kpi){kpj) / (pipj)) for all 
j I. The above factorization is directly related to the non-Abehan generalization of 
the Gribov theorem in Eq. (52). 

The form in which we present Eq. (54) corresponds to a covariant gauge for the 
gluon with momentum k. In this region for j ^ I we have PjPi/kpj ~ Ei/uj, where Ei 
is the energy of the particle with momentum pi and uu the frequency of the emitted 
gauge boson. Using the conservation of the total non-Abelian group charge: 

n 

^r«(i)M(pi,...,p,-,...,p„;fcl) = (55) 
i=i 

we can reduce the double sum over the gauge boson insertions in Eq. (54) to a single 
sum over external legs. In addition it is convenient to use the Sudakov parametrization 
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Figure 6: Feynman diagrams contributing to the infrared evolution equation (54) for 
a process with n external legs. In a general covariant gauge the virtual gluon with 
the smallest value of k± is attached to different external lines. The inner scattering 
amplitude is assumed to be on the mass shell. 



analogously to Eq. (18) and to replace the variable u by k\ according to Eq. (26). 
The infrared evolution equation then takes on the form: 

2gl " dk\ rr^H\k±\/mi, 1) dv 



A1(pi,...,P„;/X^) = A4Born(pi,-,Pn) - 77%X^ / -2^ / 



(47r)2 ^ V J\k^\/V^ V 
X CiMipu...,Pn;kl) , (56) 

where Ci is the eigenvalue of the Casimir operator T°'{l)T"'{l) {Ci = Ca for gauge bosons 
in the adjoint representation of the gauge group SU{N) and Ci = Cp for fermions in 
the fundamental representation). 

The differential form of the infrared evolution equation follows immediately from 
(56): 

= X(/x )M(pi,...,p„;// ), (57) 

where 
with 

As in the Abelian case, Wi is the probability to emit a soft and almost collinear 
gauge boson from the particle / with mass mi, subject to the infrared cut-off on 
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the transverse momentum. Note again that the cut-off is not taken to zero. The 
function r is determined by (28) for arbitrary values of the ratio mi/ ji. To logarithmic 
accuracy, we obtain from (59): 



91og(/x2) 




Qlog 



max(/x^, mf) 



s 



(60) 



The infrared evolution equation (57) should be solved with an appropriate initial con- 
dition. In the case of large scattering angles, if we choose the cut-off to be the large 
scale s then clearly there are no Sudakov corrections. The initial condition is therefore 



and the solution of (57) is thus given by the product of the Born amplitude and the 
Sudakov form factors: 



Therefore we obtain an exactly analogous Sudakov exponentiation for the gauge group 
SU{N) to that for the Abelian case, see (30). Theories with semi-simple gauge groups 
can be considered in a similar way. 

2.4 Subleading corrections from splitting functions 

At high energies, where particle masses can be neglected, the form of soft and 
coUinear divergences is universal. In this regime it is then appropriate to employ 
the formalism of the Altarelli-Parisi approach [59] and to calculate the corresponding 
splitting functions. We will do so below only for the virtual case. An important 
observation in this connection is that at high energies, all subleading terms are either 
of the coUinear or the RG type. This can be seen as follows: 

The types of soft, i.e. \k^\ <S ^/s, divergences in loop corrections with massless 
particles, unlike the coUinear logarithms, can be obtained by setting all k dependent 
terms in the numerator of tensor integrals to zero (since the terms left are of the 
order of the hard scale s). Thus it is clear that the tensor structure which emerges 
is that of the inner scattering amplitude in Fig. 8 taken on the mass-shell, times a 
scalar function of the given loop correction. In the Feynman gauge, for instance, we 
find for the well known vertex corrections the familiar three-point function Cq and for 
higher point functions we note that in the considered case all infrared divergent scalar 
integrals reduce to Co multiplied by factors of j etc.. The only infrared divergent three 
point function is given by 



M{pi,...,Pn;s 



) = MBovniPl, ■■■,Pn), 



(61) 



M{pi,...,p„;il'^) = A^Born(Pl,-,Pn)exp --^W,(s,/X^) 



(62) 





(63) 



kl>n^ (27r)4 (A;2 + ie){k'^ + 2pjk + ie){k^ - 2pik + ie) 
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The function Cq{s/ij?') is fastly converging for large k\ and we are interested here in 
the region ji^ ^ s in order to obtain large logarithms. Then logarithmic corrections 
come from the region fc^ ^ s|f | ^ s (the strong inequalities give DL, the simple 
inequalities single ones) and we can write to logarithmic accuracy: 



poo poo poo 



1 



{suv — k']_ + ie){suv — k\+ su + ie){suv — k\ — sv + is) 

Sm'^ du dv ,,2/,/, 2 2^^/ ,2n 

-7-— r / — / — / dk^eik^ - iJi^)5isuv -k\) 
2(27r)4 J-i su J-i sv J-oc ^ ^ ^ ^ ^ ^ 

i du dv ^, 2\ 

2(4^y-i-IiV^^^""-^) 

i du dv ^, 9, 

/ — / — 9(suv - IjL^) 

Jo u Jo V 



(47r)2s 
i 

(4:7lYs JiH U Jt!- V 

\ / s su 



1 du dv 



I , 2 s 



log' ^ (64) 



2(47r)2s i_i 

Thus, no single soft logarithmic corrections arc present in Co(s/yU^). In order to see 
that this result is not just a consequence of our regulator, we repeat the calculation 
for a fictitious gluon mass^. In this case we have 

Co(s/A') = j (27r)4 (F - A2 + te){k^ + 2pjk + te){k^ - 2pik + te) 

It is clear that Co(s/A^) contains soft and coUinear divergences (A- || pj^i) and is regu- 
lated with the cutoff A, which plays the role of /j, in this case. Integrating over Feynman 
parameters we find: 

We are only interested here in the real part of loop corrections of scattering amplitudes 
since they are multiplied by the Born amplitude and the imaginary pieces contribute to 
cross sections at the next to next to leading level as mentioned above. In fact, the minus 
sign inside the double logarithm corresponds precisely to the omitted principle value 



^Note that this regulator spoils gauge invariance and leads to possible inconsistencies at higher 
orders. Great care must be taken for instance when a three gluon vertex is regulated inside a loop 
integral. 
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Figure 7: In an axial gauge, all collinear logarithms come from corrections to a partic- 
ular external line (depending on the choice of the four vector satisfying n^AI, = 0) 
as illustrated in the figure. In a covariant gauge, the sum over all possible insertions 
is reduced to a sum over all n-external legs due to Ward identities. Overall, these 
corrections factorize with respect to the Born amplitude. 

contribution of Eq. (25) in the previous calculation. Thus, no single soft logarithmic 
correction is present in the case when particle masses can be neglected. 

This feature prevails to higher orders as well since it has been shown that also in 
non-Abehan gauge theories the one- loop Sudakov form factor exponentiates [50]- [58]. 

In case we would keep mass-terms, even two point functions, which in our scheme 
can only yield collinear logarithms, would contain a soft logarithm due to the mass- 
renormalization which introduces a derivative contribution (see for instance Ref. [60]). 
In conclusion, all leading soft corrections are contained in double logarithms (soft 
and collinear), and subleading logarithmic corrections in a massless theory, with all 
invariants large (s^y = '^PjPi ~ C'(s)) compared to the infrared cutoff, are of the 
collinear type or renormalization group logarithms. 

The universal nature of collinear type logarithmic corrections can then easily be 
seen in an axial gauge where collinear logarithms are related to corrections on a partic- 
ular external leg depending on the choice of the four vector rii, [61, 62, 63, 64]. A typical 
diagram is depicted in Fig. 7. In a general covariant gauge this corresponds (using 
Ward identities) to a sum over insertions in all n external legs [28]. We can therefore 
adopt the strategy to extract the gauge invariant contribution from the external line 
corrections on the invariant matrix element at the subleading level. The results of 
the above discussion are thus important in that they allow the use of the AltareUi- 
Parisi approach to calculate the subleading contribution to the evolution kernel of Eq. 
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(57). We are here only concerned with virtual corrections and use the universality of 
the splitting functions to calculate the subleading terms. For brevity wc discuss both 
scalar as well as conventional QCD simultaneously. In each case one only needs to 
switch off the other type of field in the fundamental representation to obtain the case 
of interest. The /3-function in both cases differs in the non-glue part but since this 
difference is of no consequence in our later discussion we don't distinguish between the 
two. For the purpose of calculating SL virtual corrections we use the virtual quark, 
scalar quark and gluon contributions to the splitting functions P^^{z), Pj^ and Pjg{z) 
describing the probability to emit a soft and/or coUinear virtual particle with energy 
fraction z of the original external line four momentum. The infinite momentum frame 
corresponds to the Sudakov parametrization with lightlike vectors. In general, the 
splitting functions Pba describe the probability of finding a particle B inside a particle 
A with fraction z of the longitudinal momentum of A with probability Vba to first 
order [59]: 

dVBA{z) = -^PsAdt (67) 

ZTT 

where the variable t — log ^ for our purposes. It then follows [59] that 



|2 



JT) ^^ Ois z{l - z)^^\Va^b+c\- ,2 /«oN 
dVBA{z)^— 2^ -2 '-dlogk^ (68) 

spins -L 



where Va >b+c denotes the elementary vertices and 
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Z{1 - z) ^ \Va^B+C\ .^qx 

Pba{z) = 2_. p (69) 

spins -L 

The upper bound on the integral over dk\ in Eq. (68) is s and it is thus directly related 
to dt. Regulating the virtual infrared divergences with the transverse momentum cutoff 
as described above, we find the virtual contributions to the splitting functions for 
external quark, scalar quarks and gluon lines: 

= Cf (^-2 log ^ + 3^5(1 -z) (70) 

PZiz)^CF(^-2\og^+ 4^6(1 -z) (71) 

= Ca (-2 log ^ + ^/So'^'^^) '^(l - z) (72) 

The functions can be calculated directly from loop corrections to the elementary pro- 
cesses [65, 66, 67, 38] and the logarithmic term corresponds to the leading kernel of 
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section 2.3. We introduce virtual distribution functions which include only the effects 
of loop computations. These fulfill the AltareUi-Parisi equations^ 



dq{z,t) 

dt 
ds{z, t) 

dt 
dg{z,t) 

dt 



27r Jz y 



■g(z/y,t)PM 



(73) 
(74) 
(75) 



The splitting functions are related by Pba = Pba + Pba-> where R denotes the con- 
tribution from real gauge boson emission^. Pba is free of logarithmic corrections and 
positive definite. The subleading term in Eq. (72) indicates that the only subleading 
corrections in the pure glue sector are related to a shift in the scale of the couphng. 
These corrections enter with a different sign compared to the conventional running 
coupling effects. For fermion and scalar external lines there is an additional subleading 
correction from coUinear terms which is not related to a change in the scale of the 
coupling. 

Inserting the virtual probabilities of Eqs. (70), (71) and (72) into the Eqs. (73), 
(74) and (75) we find: 



g(l,t) 
s{l,t) 



?oexp 
So exp 
5(0 exp 



q;.,C, 



— (^log ^-31og-^^ 
[log J.-^^og-^ 

UsCa a 2 S 4 QCD, S 



1^ 



(76) 
(77) 
(78) 



where (3'^^^ = %Ca - iTpUf with Ca = 3, Cf = 4/3 and Tp = ^ 
describe the total contribution for the emission of virtual particles 
all invariants large compared to the cutoff /x, to the densities q{z, t), s{z, t) and g(z, t). 
The normalization is on the level of the cross section. For the invariant matrix element 
we thus find at the subleading level for processes with n external lines: 

Miih Pn-g,. fi) = M{pi, ...,Pn, gs) x 



These functions 
i.e. z — 1), with 



*Note that the off diagonal splitting functions Pqg and Pgg etc. do not contribute to the virtual 
probabilities to the order we are working here. In fact, for virtual corrections there is no need to 
introduce off-diagonal terms as the corrections factorize with respect to the Born amplitude. The 
normalization of the Eqs. (70), (71) and (72) corresponds to calculations in two to two processes on 
the cross section level with the gluon symmetry factor ^ included. The results, properly normalized, 
arc process independent. 

^Pqq was first calculated by V.N. Gribov and L.N. Lipatov in the context of QED [68, 69]. 



25 



Figure 8: Feynman diagrams contributing to the infrared evolution equation (84) for 
a process with n external legs. In a general covariant gauge the virtual gluon with 
the smallest value of fcj^ is attached to different external lines. The inner scattering 
amplitude is assumed to be on the mass shell. 



1 -X^. , o. 1 . 1 " 
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exp I -- E - ^ E "^s^s, - ^ E I (79) 

^ 3=\ ^ i=l ^ l=\ 



with ^5 + ^3 + ^5 = n, and 



»'A/'^) = ^(^V^-31og^j (80) 
W,(^./<^)=^^(W^-41ogij) (81) 
Z-^) = ^ f V 4. - ^,3?=° log ^\ (82) 



The functions VF^, Wg and VFg correspond to the probability of emitting a virtual soft 
and/or coUinear gauge boson from the particle g, g subject to the infrared cutoff \x. 
Typical diagrams contributing to Eq. (79) in a covariant gauge are depicted in Fig. 8. 
In massless QCD there is no need for the label Wq^^ Wg^ or Wg^, however, we write it 
for later convenience. The universality of the splitting functions is crucial in obtaining 
the above result. 



2.5 Anomalous scaling violations 

The solution presented in Eq. (79) determines the evolution of the virtual scat- 
tering amplitude ^A{pl, ...,Pn, gs, /J') for large energies at fixed angles and subject to 
the infrared regulator fi. In the massless case there is a one to one correspondence 
between the high energy limit and the infrared limit as only the ratio s/jj,^ enters as a 
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dimensionless variable [70, 71]. Thus, we can generalize the Altarelli-Parisi equations 

(73), (74) and (75) to the invariant matrix element in the language of the renormaliza- 
tion group. For this purpose, we define the infrared singular (logarithmic) anomalous 
dimensions 

TM . ; TAt) . ; T.M . (83) 

Infrared divergent anomalous dimensions have been derived in the context of renor- 
malization properties of gauge invariant Wilson loop functionals [72]. In this context 
they are related to undifferentiable cusps of the path integration and the cusp angle 
PjPi/lJ>^ gives rise to the logarithmic nature of the anomalous dimension. In case we use 
off-shell amplitudes, one also has contributions from end points of the integration [72]. 
The leading terms in the equation below have also been discussed in Refs. [73], [74] 
and [75] in the context of QCD. With these notations we find that Eq. (79) satisfies 

+ns (r,{t) + ^7,)) Mipi, ...,pn,gs,l^) = (84) 

to the order wc arc working here and where A4{pi, ...,Pn, Qs, fi) is taken on the mass- 
shell. The difi'crcnce in the sign of the derivative term compared to Eq. (57) is 
due to the fact that instead of differentiating with respect to log//^ we use logs//x^. 
The anomalous dimensions are given by 7p = — ~ (ll^^ ~ ^^-p)' 7* ~ 

—Cp^ and 7g = — Cf|-. As mentioned above in pure scalar QCD the /3-function 
differs in the non-glue part from f3'^'~"^. The quark- antiquark operator anomalous 
dimension 7^ or in scalar QCD 7^ enter even for massless theories as the quark antiquark 
operator leads to scaling violations through loop effects since the quark masslessness 
is not protected by gauge invariance and a dimensionful infrared cutoff needs to be 
introduced. Thus, although the Lagrangian contains no mipip or m^cffcf) term, quantum 
corrections lead to the anomalous scaling violations in the form of 7^ or 7^. The factor 
I occurs since we write Eq. (84) in terms of each external line separately^. For the 
gluon, the scaling violations due to the infrared cutoff are manifest in terms of an 
anomalous dimension proportional to the /3-function since the gluon mass is protected 
by gauge invariance from loop corrections. Thus, in the bosonic sector the subleading 
terms correspond effectively to a scale change of the coupling. Fig. 9 illustrates the 
corrections to the external quark- antiquark lines from loop effects. 

Except for the infrared singular anomalous dimension (Eq. (83)), all other terms 
in Eq. (84) are the standard contributions to the renormalization group equation for 

^In case of a massive theory, we could, for instance avoid the anomalous dimension term 7^ by 
adopting the pole mass definition. In this case, however, we would obtain terms in the wave function 
renormalization, and in any case, the one to one correspondence between UV and IR scaling, crucial 
for the validity of Eq. (84), is violated. 
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Figure 9: The two counterterms contributing to the quark anomalous dimension 7^ = 
aiogTi" (-^99 + ^2) or 7s = giog^2 (Ss's + ^2)- Here JI denotes the MS dimensional 
regularization mass parameter. Due to divergences in loop corrections there are scaling 
violations also in the massless theory. 

S-matrix elements [76]. In QCD, observables with infrared singular anomalous dimen- 
sions, regulated with a fictitious gluon mass, are ill defined due to the masslessness 
of gluons. In the electroweak theory, however, we can legitimately investigate only 
virtual corrections since the gauge bosons will require a mass. Eq. (84) will thus be 
very useful in section 3. 

2.6 Renormalization group corrections 

In this section we review the case of higher order RG-corrections in unbroken 
gauge theories like QCD following Ref. [77]. Explicit comparisons with higher or- 
der calculations for the on-shell Sudakov form factor revealed that the relevant RG 
scale in the respective diagrams is indeed the perpendicular Sudakov component 
[78, 79, 80, 81, 82]. It should be noted, however, that in particular for the mas- 
sive cases, the discussion below is not vahd close to any of the thresholds. In such 
cases it is useful to consider "physical renormalization schemes" such as discussed in 
Ref. 3,mV2,mV3itemV,mV2,mV3, which display a gauge invariant, continuous and 
smooth flavor threshold behavior with automatic decoupling of heavy particles. For 
our purposes here, we give correction factors for each external line below. The univer- 
sal nature of the higher order SL-RG corrections can be seen as follows. Consider the 
gauge invariant fermionic part (~ ri/) as indicative of the full P^'^^ term (replacing 

% = ^ (tI^^ - Po^^))- In order to lead to subleading, i.e. O (al^log^"-^ this 
loop correction must be folded with the exchange of a gauge boson between two exter- 
nal lines (producing a DL type contribution) like the one depicted in Fig. 10. Using 
the conservation of the total non-Abelian group charge, i.e. Eq. (55), the double sum 
over all external insertions j and I is reduced to a single sum over all n external legs. 
Thus these types of corrections can be identified with external lines at higher orders. 
The same conclusion is reproduced by the explicit pole structure of MS renormalized 
scattering amplitudes at the two loop level in QCD [86]. The results presented in 
Ref. [86] have been confirmed recently by exphcit massless two loop QCD calculations 
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Figure 10: A QED diagram at the two loop level yielding a SL-RG correction. The 
explicit result obtained in Ref. [87] for the case of equal masses relative to the Born 
amphtude was — 3^ = i^/^o if;;^ ^og This result is reproduced exactly 

by including a running coupling into the one loop vertex correction diagram. The 
argument of the coupling must depend on the component of the loop momentum 
(going into the fermion loop) which is perpendicular to the external fermion momenta. 
In QCD, although more diagrams contribute, the net effect is just to replace — > 
Pq'^^ in the above expression. 

[98, 99, 100, 101, 102, 103]. In addition, from the expression in Ref. [86] it can be 
seen that the SL-RG corrections are independent of the spin, i.e. for both quarks and 
gluons the same running coupling argument is to be used. This is a consequence of 
the fact that these corrections appear only in loops which can yield DL corrections on 
the lower order level and as such, the available DL phase space is identical up to group 
theory factors. We begin with the virtual case. 

2.6.1 Virtual corrections 

The case of virtual SL-RG corrections for both massless and massive partons has 
been discussed in Ref. [88] with a different Sudakov paramctrization. Below we show 
the identity of both approaches. The form of the corrections is given in terms of the 
probabilities Wi^ (s, /i^). To logarithmic accuracy, they correspond to the probability 
to emit a soft and/or coUinear virtual parton from particle i at high energies subject 
to an infrared cutoff /i. At the amplitude level all expressions below are universal for 
each external line and exponentiate according to 

M{pi,...,Pn,gs,l^) ^ MBorniPl, ■■■,Pn, 9s) (^—J^Wi^ (85) 

where n denotes the number of external lines. We begin with the massless case. 
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Massless QCD 

In the following we denote the running Q CD-coupling by 

<kl) = ""-^ ^ ^ ^Al^ (86) 

l + ^/?o^"°log^ 1 + clog^ 

Up to two loops the massless /3-function is independent of the chosen renormalization 
scheme and is gauge invariant in minimally subtracted schemes to all orders [89]. These 
features will also hold for the derived renormalization group correction factors below 
in the high energy regime. The scale ji denotes the infrared cutoff on the exchanged 
fex between the external momenta where the Sudakov decomposition is given by 

k = vpi + upj + k±, such that Pjk± — pik± = 0. The cutoff /i serves as a a lower limit 
on the exchanged Euclidean component fc^ = — A;^ > as in the previous sections. 
In order to avoid the Landau pole we must choose /i > Aqcd- Thus, the expressions 
given in this section correspond for quarks to the case where m <^ /i. For arbitrary 
external lines we then have 



The RG correction is then described by including the effect of the running coupling 
from the scale //^ to s according to [78, 79, 80, 81, 82] (see also discussions in Refs. 
[88, 90]): 

Ci dkl /•! dv n,(//-) 




^'-/^ ^ 1 + c log 



1 , si, aJii'^) \ 1 , a^iiu'^) 1 . s 

c'°vrw-v^?'°^ii/l '''' 

where Cj = Ca for gluons and Ci = Cp for quarks. For completeness we also give 
the subleading terms of the external line correction which is of course also important 
for phenomenological applications. The terms depend on the external line and the 
complete result to logarithmic accuracy is given by: 

„o = ^M)£a 1 1 ,„g 4 _ i) + i log ^ 

y ' 27r [c ii- \ OLs\s) ) & OLs\s) 

-7^/?o^""log4^ (89) 
Ca J 
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It should be noted that the subleading term in Eq. (89) proportional to Pq is 

not a conventional renormaUzation group corrections but rather an anomalous scaling 
dimension, and enters with the opposite sign [37] compared to the conventional RG 
contribution (see section 2.5). 



Massive QCD 

Here we give results for the case when the infrared cutoff /i. -C m, where m denotes 
the external quark mass. We begin with the case of equal external and internal line 
masses: 



Equal masses 



Following Ref . [88] , we use the gluon on-shell condition suv = k\ to calculate the 
integrals. We begin with the correction factor for each external massive quark line. 
Following the diagram in Fig. 2 we find: 
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The //-dependent terms cancel out of any physical cross section (as they must) when 
real soft Bremsstrahlung contributions are added and c = as{m^)j3Q'^^ /n for massive 
quarks. In order to demonstrate that the result in Eq. (91) exponentiates, we calcu- 
lated in Ref. [88] the explicit two loop renormalization group improved massive virtual 
Sudakov corrections, containing a different "running scale" in each loop. It is of course 
also possible to use the scale k\ directly. In this case we have according to the diagram 
in Fig. 3: 
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which is the identical result as in Eq. (91). For completeness we also give the subleading 
terms of the pure one loop form factor which is again important for phenomenological 
applications. The complete result to logarithmic accuracy is thus given by: 
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(93) 



2 n 

For m — ji Eq. (93) agrees with Eq. (90) in the previous section for massless quarks. 



Unequal masses 

In this section we denote the external mass as before by m and the internal mass 
by nii and thus, the constant c = Q;s(mf )/3^'"°/7r. We consider only the case at high 
energies taking the first two families of quarks as massless. The running of all light 
flavors is implicit in the Uf term of the function. The result is then given by: 
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(94) 



It is evident that the effect of unequal masses is large only for a large mass splitting. 
In QCD, we always assume scales larger than Aqcd and with our assumptions we have 
only the ratio of rrit/mb leading to significant corrections. 
The full subleading expression is accordingly given by: 
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For m = rrii Eq. (95) agrees with Eq. (93) in the previous section for equal mass 
quarks. 
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If we want to apply the above result for the case of QED corrections later, then 
there is no Landau pole (at low energies) and we can have large corrections of the form 
mfe/me etc. In this case the running coupling term is given by 

e\kl) = ^ (96) 

and instead of Eq. (95) we have: 



i log ^^^^ 1 - i ^ E OjiVJ log - I log 4 - log (97) 

and where c = E%i Q^N^.. 

2.6.2 Real gluon emission 

We discuss the massless and massive case separately since the structure of the diver- 
gences is different in each case. For massive quarks we discuss two types of restrictions 
on the experimental requirements, one in analogy to the soft gluon approximation. 
The expressions below exponentiate on the level of the cross section, i.e. for observable 
scattering cross sections they are of the form 

d(7(pi, ...,p„,g's,//expt) = daBoTn{,Pl,-;Pn,gs) X 

exp [W,,n {s, ix\ i^l^,) - W^y {s, ^?)] | (98) 

where the sum in the exponential is independent of ^ and only depends on the cutoff 
A*expt defining the experimental cross section. We begin with the massless case. 

Emission from massless partons 

In this section we consider the emission of real gluons with a cutoff k± < /iexpt, 
related to the experimental requirements. For massless partons we have at the DL 
level: 
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and thus for the RG-improved correction: 

Cj r^oxpt dk^ duo as{fi'^ 



This expression depends on /i as it must in order to cancel the infrared divergent virtual 
corrections. In fact the sum of real plus virtual corrections on the level of the cross 
section is given by 
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and thus independent of fi. The full expressions to subleading accuracy are thus: 
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All divergent (//-dependent) terms cancel when the full virtual corrections are added. 



Emission from massive quarks 

In the case of a massive quark, i.e. ji <^ m, the overall infrared divergence is 
not as severe. This means we can discuss different requirements which all have the 
correct divergent pole structure canceling the corresponding terms from the virtual 
contributions. We divide the discussion in two parts as above. 
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Equal masses 

The constant c = Q;s(m^)/3^*^°/7r below. We have the following expression without 
a running coupling: 
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m 



If we want to employ a restriction analogously to the soft gluon approximation, we 
find independently of the quark mass [37, 30] : 
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In all cases above we have not taken into account all subleading coUinear logarithms 
related to real gluon emission. In order to now proceed with the inclusion of the 
running coupling terms it is convenient to first consider only the DL phase space in 
each case. Thus we find 
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The full subleading expressions are thus given by 
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and 
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In case we also impose a cut on the integration over u we have independently of the 
relation between m and Hexpt assuming only <^ s: 
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This expression agrees with the result obtained in Ref. [88] where the gluon on-shell 
condition fc^ = suv was used and one integral over one Sudakov parameter was done 
numerically. In Ref. [88] it was also shown that the RG-improved virtual plus soft 
form factor also exponentiates by explicitly calculating the two loop RG correction 
with each loop containing a running coupling of the corresponding k\. 

The full subleading expression for the RG-improved soft gluon emission correction 
is thus given by 
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for the equal mass case. The case of different external and internal masses is again 
important for applications in QED and will be discussed next. 



Unequal masses 

While the gluonic part of the /3-function remains unchanged we integrate again 
only from the scale of the massive fermion which is assumed to be in the perturbative 
regime. For applications to QED, however, we need the full expressions below. Here 
we discuss only the case analogous to the soft gluon approximation. Considering again 
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only the high energy scenario we have for the case of an external mass m and a fermion 
loop mass mf 
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This expression agrees with the result obtained in Eq. (110) for the case rrii = m. 

The full subleading expression for the RG-improved soft gluon emission correction 
is thus given by 
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As mentioned above, this expression is more useful for applications in QED or if the 
mass ratios arc very large. In QED we have again the running coupling of the form 
given in Eq. (96), and Eq. (113) becomes 
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where again c = — Sj^i Q'j^c- 'This concludes the discussion of SL-RG effects in 
QCD. As a side remark we mention that for scalar quarks, the same function appears 
as for fcrmions since the DL-phase space for both cases is identical. Only (3q differs in 
each case. 
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3 Broken gauge theories 

In the following we will apply the results obtained in the previous sections to 
the case of spontaneously broken gauge theories. It will be necessary to distinguish 
between transverse and longitudinal degrees of freedom. The physical motivation in 
this approach is that for very large energies, s S> = M^, the electroweak theory is 
in the unbroken phase, with effectively an SU{2) x f/(l) gauge symmetry as described 
by the high energy symmetric part of the Lagrangian in Eq. (11). We will calculate 
the corrections to this theory and use the high energy solution as a matching condition 
for the regime for values of /i < M. 

We begin by considering some simple kinematic arguments for massive vector 
bosons. A vector boson at rest has momentum k'^ = (M, 0, 0, 0) and a polarization 
vector that is a linear combination of the three orthogonal unit vectors 

ei = (0,l,0,0) , ei = (0,0,l,0) , eg = (0, 0, 0, 1) . (115) 

After boosting this particle along the 3-axis, its momentum will be k'^ = {Ek, 0, 0, k). 
The three possible polarization vectors are now still satisfying: 

k^e'< = , £| = -1 . (116) 

Two of these vectors correspond to ei and 62 and describe the transverse polarizations. 
The third vector satisfying (116) is the longitudinal polarization vector 

elik)^ik/M,0,0,Ek/M) (117) 

i.e. s1{k) = k^ /M + 0{M/Ek) for large energies. These considerations illustrate that 
the transversely polarized degrees of freedom at high energies are related to the massless 
theory, while the longitudinal degrees of freedom need to be considered separately. 

Another manifestation of the different high energy nature of the two polarization 
states is contained in the Goldstone boson equivalence theorem. It states that the 
unphysical Goldstone boson that is "eaten up" by a massive gauge boson still controls 
its high energy asymptotics. A more precise formulation is given below in section 3.3. 

Thus we can legitimately use the results obtained in the massless non-Abelian 
theory for transverse degrees of freedom at high energies and for longitudinal gauge 
bosons by employing the Goldstone boson equivalence theorem. 

Another difference to the situation in an unbroken non-Abelian theory is the mixing 
of the physical fields with the fields in the unbroken phase. These comphcations are 
especially relevant for the Z-boson and the photon. 

3.1 Fermions and Transverse degrees of freedom 

The results we obtain in this section are generally valid for spontaneously broken 
gauge theories, however, for definiteness we discuss only the electroweak Standard 
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Model. The physical gauge bosons are thus a massless photon (described by the field 
Ai,) and massive and Z bosons (described correspondingly by fields and Z^).: 

Wt^^i^l^iWl) (118) 

Z^ = cos d^Wl + sin Q^B^ (119) 
A^ ^- sin d^Wl + cos d^B^ (120) 

Thus, amplitudes containing physical fields will correspond to a linear combination of 
the massless fields in the unbroken phase. The situation is illustrated schematically for 
a single gauge boson external leg in Fig. 11. In case of the bosons, the corrections 
factorize with respect to the physical amplitude. 

In the general case let us denote physical particles (fields) by / and particles 
(fields) of the unbroken theory by u. Let the connection between them be denoted by 
/ = Yl,u C^'^u, where the sum is performed over appropriate particles (fields) of the 
unbroken theory as in Eqs. (118), (119) and (120). Note that, in general, physical 
particles, having definite masses, don't belong to irreducible representations of the 
symmetry group of the unbroken theory (for example, the photon and Z bosons have 
no definite isospin). On the other hand, particles of the unbroken theory, belonging 
to irreducible representations of the gauge group, have no definite masses. Then for 
the amplitude A1-^i' '-^"(j9i, ...,Pn', ^'^) with n physical particles fi with momenta pj and 
infrared cut-off /x^, the general case for virtual corrections is given by 

n 

At^^'-^"(pi,...,p„;/x^)= n^^'"^-^"''-""bi--Pn,/x') (121) 

U1,...U„ j=l 

On the one loop level and to subleading accuracy, Eq. (121) must also include the 
correct counterterms for the commonly chosen on-shell scheme. In this scheme the on- 
shell photon, for instance, does not mix with the Z-boson, thus including all mixing 
effects into the massive neutral Z-boson sector. More details are given below. In the 
following we give results only for the amplitudes 7W"i' - ""(]9i, ...,Pn, At^) keeping in mind 
that in general the physical amplitudes must be obtained via Eq. (121). For fermions, 
transverse W^, longitudinal gauge bosons or Higgs bosons, no linear combination 
arises, i.e. the universal corrections below factorize automatically with respect to the 
physical Born process. Only photons and Z-bosons are affected by this complication 
for the obvious reasons discussed above. To logarithmic accuracy, all masses can be 
set equal: 

Mz ~ Mw ~ ^Higgs ~ M 

and the energy is considered to be much larger, ^/s ^ M. The left and right handed 
fermions are correspondingly doublets (T = 1/2) and singlets (T = 0) of the SU{2) 
weak isospin group and have hypercharge Y related to the electric charge Q, measured 
in units of the proton charge, by the Gell-Mann-Nishijima formula Q = + Y /2. 
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Figure 11: The schematic corrections to external gauge boson emissions in terms of 
the fields in the unbroken phase of the electroweak theory. There are no mixing terms 
between the Wj^ and Bi, fields for massless fermions. We denote cos 9^ by and 
sin 9^ by Sw For final states, the corrections factorize with respect to the physical 
amplitude. In general, one has to sum over all fields of the unbroken theory with each 
amplitude being multiplied by the respective mixing coefficient as given in Eq. (121). 
At one loop, also the renormalization conditions must be included. 
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The value for the infrared cutoff fj, can be chosen in two different regimes (see Fig. 
1): I ^/s fi ^ M and II ^ M. The second case is universal in the sense that it 
does not depend on the details of the electroweak theory and will be discussed below. 
In the first region we can neglect spontaneous symmetry breaking effects (in particular 
terms connected to the v.e.v.) and consider the theory with fields B,^ and as given 
by >Csymm in Eq. (11). One could of course also calculate everything in terms of the 
physical fields, however, we emphasize again that in this case we need to consider the 
photon also in region I). The omission of the photon would lead to the violation of 
gauge invariance since the photon contains a mixture of the B^, and Wj^ fields. 

In region I), the renormalization group equation (or generalized infrared evolution 
equation) (84) in the case of all rrii < M reads^ 

xM^{pi,...,Pn,g,g',iJ>) = (122) 

where the index _L indicates that we consider only Ug transversely polarized external 
gauge bosons with nw + = Ug and Uf denotes the number of external fermion lines. 
The two /9-functions are given by: 

PiMp')) =|f^^l»-A^ (123) 

with the one-loop terms given by: 

Po = - ^rigen - ^rih , P'o = —rigen - (125) 

where n^en denotes the number of fermion generations [91, 92] and Uh the number of 
Higgs doublets. Eq. (122) describes the one loop RG corrections correctly. At higher 
orders the subleading RG terms must be included according to the discussion in section 
3.5. The infrared singular anomalous dimensions read 

= f^^^m + 1) + ^ f^yi t (126) 



47r ' ' 47r 



^Note, that the ampUtude on the right hand side is in general a hnear combination of fields in 
the unbroken phase according to Eq. (121). In addition, in the electroweak theory matching will be 
required at the scale M and often on-shell renormalization of the couplings e and sin 9^ is used. In 
this case one has additional complications in the running coupling terms due to the different mass 
scales involved below M. Details are presented in section 3.5. 
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where 7] and 1^ are the total weak isospin and hypercharge respectively of the particle 
emitting the soft and coUinear gauge boson. Analogously, 
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where the last two terms only contribute for quarks of the third generation, i' denotes 
the isospin partner of i. The presence of Yukawa terms and also the Higgs contribution 
to the /3-functions in Eq. (125) are remnants of the spontaneous broken symmetry 
which leads to differences even in the transverse sector compared to unbroken gauge 
theories as is obvious from the form of £symm in Eq. (11). In terms of the corresponding 
logarithmic probabilities we thus have the following expression for fermions from the 
virtual splitting function approach: 
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For external transversely polarized gauge bosons: 
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The initial condition for Eq. (122) is given by the requirement that for the infrared 
cutoff fj,^ — s we obtain the Born amplitude. The solution of (122) is thus given by 

M^{pi, ...,Pn, g, 9, IJ) = A^Born(Pl, -^Pn, 9, 9) 
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where we neglect RG corrections for now. These will be discussed thoroughly in section 
3.5. nw and ne denote the number of external W and B fields respectively. The 
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SU (2) xU{l) group factors in the exponential can be written in terms of the parameters 
of the broken theory as follows: 



g%{T, + l)+g'' 
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cos^ 9n 



where the three terms on the r.h.s. correspond to the contributions of the soft photon 
(interacting with the electric charge Cj = QigsvuO^)^ the and the Z bosons, re- 
spectively. Although we may rewrite solution (130) in terms of the parameters of the 
broken theory in the form of a product of three exponents corresponding to the ex- 
changes of photons, and Z bosons, it would be wrong to identify the contributions 
of the diagrams without virtual photons with this expression for the particular case 
e| = 0. This becomes evident when we note that if we were to omit photon lines then 
the result would depend on the choice of gauge, and therefore be unphysical. Only for 
6y, = 0, where the photon coincides with the B gauge boson, would the identification 
of the e? term with the contribution of the diagrams with photons be correct. 

We now need to discuss the solution in the general case. In region I) we calculated 
the scattering amplitude for the theory in the unbroken phase in the massless limit. 
Choosing the cutoff ^ in region II), ^ ^ M, we have to only consider the QED 
contribution. In this region we cannot necessarily neglect all mass terms, so we need 
to discuss the subleading terms for QED with mass effects. If <^ ji, the results 
from massless QCD can be used directly by using the Abehan hmit Ci? = 1. In case 
-C mi we must use the well known next to leading order QED results, e.g. [93], and 
the virtual probabilities take the following form for fermions: 
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Note, that in the last equation the full subleading coUinear logarithmic term [60] is 
used in distinction to Ref. [93]. In the explicit two loop calculation presented in Ref. 
[94] it can be seen that the full collinear term also exponentiates at the subleading 
level in massive QED. For bosons we have analogously: 
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In addition we have collinear terms for external on-shell photon lines from fermions 
with mass rrij and electromagnetic charge ej up to scale M: 
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Note that automatically, w^X^'^^ -^^) — 0- one loop order, this contribution cancels 

against terms from the rcnormalization of the QED coupling up to scale M . For 
external Z-bosons, however, there are no such collinear terms since the mass is large 
compared to the rrii. Thus, the corresponding RG- logarithms up to scale M remain 
uncanceled. 

The appropriate initial condition is given by Eq. (130) evaluated at the matching 
point ji — M. Thus we find for the general solution in region II): 
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The last equality holds for ^ mj <^ M and we have replaced the matching scale M by 
rrit in the Yukawa enhenced subleading terms since the coefficients are unambiguously 



44 



determined and the argument in the corresponding logarithm must be rrit [38, 29]. It 
is important to note again that, unhke the situation in QCD, in the electroweak theory 
we have in general different mass scales determining the running of the couplings of the 
physical on-shell renormalization scheme quantities. We have written the above result 
in such a way that it holds for arbitrary chiral fermions and transversely polarized gauge 
bosons. In order to include physical external photon states in the on-shell scheme, the 
renormalization condition is given by the requirement that the physical photon does 
not mix with the Z-boson. This leads to the condition that the Weinberg rotations in 
Fig. 11 at one loop receive no RG- correct ions. Thus, above the scale M the subleading 
coUinear and RG-corrections cancel for physical photon and Z-boson states. 

Since the Yukawa enhanced terms are novel features in broken gauge theories as 
compared to the situation in QCD we use the non-Abelian generalization of the Gri- 
bov theorem in the following to prove the correctness of our splitting function ap- 
proach for specific processes. Since we are interested here in corrections to order 
O (a'^ log^"~^ , each additional loop correction to the universal subleading terms in 
the previous section must yield two logarithms, i.e. we are considering DL-corrections 
to the basic process like the inner fermion loop in Fig. 13. It is of particular importance 
that all additional gauge bosons must couple to external legs, since otherwise only a 
subleading term of order O (a^ log^"^^ jj2^ would be generated. All subleading correc- 
tions generated by the exchange of gauge bosons coupling both to external Goldstone 
bosons and inner fermion lines cancel analogously to a mechanism found in Ref. [95] 
for terms in heavy quark production in 77-collisions in a = state. Formally this 
can be understood by noting that such terms contain an infrared divergent correction. 
The sum of those terms, however, is given by the Sudakov form factor. Thus any 
additional terms encountered in intermediate steps of the calculation cancel. 

This point can be understood also from the principle of gauge invariance. At 
the two loop level for instance we have to consider diagrams of the type depicted 
in Fig. 12. They involve Vertex corrections rj,(A;^,/^, (ki — Z)^) and self energy terms 
S((A;i-/)2) with the same overall Yukawa-term structure. Writing the gauge coupling in 
the symmetric basis for clarity since we are considering a regime where s = (/ci — A;2)^ ^ 
M^, where M is the gauge boson mass. In any case, local gauge invariance is not 
violated in the SM and for heavy particles in the high energy limit, we can perform 
the calculation in a basis which is more convenient. For our purposes we need to 
investigate terms containing three large logarithms in those diagrams. Since the would- 
be Goldstone boson loops at one loop only yield a single logarithm it is clear that 
the gauge boson loop momentum I must be soft. Thus we need to show that the 
UV logarithm originating from the k integration is identical (up to the sign) in both 
diagrams. We can therefore neglect the loop momentum I inside the fermion loop. It 
is then straightforward to see that 

^j:(kl)^r^(klo,kl) (135) 
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Figure 12: Two loop Feynman diagrams yielding Yukawa enhanced logarithmic cor- 
rections to the third generation of fermions in the final sate. The Ward identity in Eq. 
(135) assures that in the Feynman gauge, the sum of all vertex and self energy diagrams 
does not lead to additional SL logarithms at the two loop level. Only corrections to 
the original one loop vertex need to be considered and lead to the exponentiation of 
Yukawa terms in the fermionic SM sector to SL accuracy. 

where the full sum of all contributing self energy and vertex diagrams must be taken. 
Thus, we have established a Ward identity for arbitrary Yukawa couplings of scalars 
to fermions and thus, the identity of the UV singular contributions. The relative sign 
is such that the generated SL logarithms of the diagrams in Fig. 12 cancel each other. 
The existence of such an identity is not surprising since it expresses the fact that also 
the Yukawa sector is gauge invariant. We are thus left with gauge boson corrections 
to the original vertices in the on-shell renormalization scheme such as depicted in Fig. 
13. At high energies we can therefore employ the non-Abelian version of Gribov's 
bremsstrahlung theorem. The soft photon corrections are included via matching as 
discussed above. 
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For the one loop process in Fig. 13, for instance, we include only corrections 
with top and bottom quarks and assume on-shell renormalization. Thus the correc- 
tions at higher orders factorize with respect to the one loop fermion amplitude and 
■A4 "Born" {pi, ■■■,Pn) = -M. Uoop{pi, Pn) ■ Notc that the latter is also independent of 
the cutoff II since the fermion mass serves as a natural regulator. In principle we can 
choose the top-quark mass to be much larger than fj, for instance. In our case we have 
for the electroweak DL corrections at the weak scale ji — M: 



-Ti(T,+ l) + 



167r2 ' ' ' ' 167r2 4 



(136) 



We now want to consider specific processes relevant at future e^e' colliders and demon- 
strate how to apply the non-Abelian version of Gribov's factorization theorem for the 
higher order corrections. The subleading corrections are then compared to the general 
splitting function result in Eqs. (128) and (130). Below we use the physical fields for 
the respective contributions. 

From the arguments of section 6.1 it is now straightforward to include also top- 
Yukawa terms for chiral quark final states. These terms occur for left handed bottom 
as well as top quark external lines. The situation for a typical Drell-Yan process is 
depicted in Fig. 13 where for the inner scattering amplitude we have two contributions. 

We neglect all terms of order (9 1^, — j . Using on-shell renormalization we find for 

the inner amplitude on the left in Fig. 13 for a right handed electron in the initial and 
a left handed bottom quark in the final state from the 0^ loop for the sum of the 7 
and Z contributions according to the Feynman rules in appendix 6.2: 

4 2 

a\ DY € Ulf / + 1 I — \ 

{him-h-hrih) DY 

{27ry {P - m), + ie){{l - hY - + ie){{l - - + ie) 

4 9 
7 P ill 

* {et\lAeK){kmh){B2,- B^,) (137) 



/ 
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The scalar functions at high energy evaluate to i?23 — = — log For the diagram 
on the right in Fig. 13 we have for the bottom again only the 0^ contribution. Here 
we find for the sum of the 7 and Z contributions: 



X 



64DY _ " -"^t^t / +1 I - 



(/L|f7^^/R) , . .DY 

(27r)" (/2 - M2 + ie){{l - - mf + ie){{l - kif - mf + ie) 



^"^tQt (e+|7,|eS)(/L|7l/R>(S23-S,^) (138) 



327r2 2sM244 
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Figure 13: Feynman diagrams yielding Yukawa enhanced logarithmic corrections to 
the third generation of fermions in the final sate. The inner scattering amplitude is 
taken on the mass shell. No DL-corrections originate from the inner loop. At higher 
orders, the sublcading corrections are given in factorized form according to the non- 
Abelian generalization of Gribov's theorem as described in the text. Corrections from 
gaTigc bosons inside the Goldstone-boson loop give only sub-sublcading contributions. 
DL-corrections at two and higher loop order are given by gauge bosons coupling to (in 
principle all) external legs as schematically indicated. 



In all cases we renormalize on-shell, i.e. by requiring that the vertex vanishes when 
the momentum transfer equals the masses of the external on-shcll lines. All on-shcll 
self energy contributions don't contribute in this scheme. For external left handed top 
quarks, the 0^ loop is mass suppressed and we only have to consider the x ^'^^ H 
corrections. They are given by replacing Qt — > 2Qt {Tf)^ and Qt — * \Qt in Eq. 
(138). It turns out that the ZxH contributions equal the corrections from the 70=*= 
and Zcf)'^ in the case of the bottom calculation. The Born amplitude is given by: 

M^l^ = t-^{Qf -Tf){eMe^){fL\Y\fR) 



sc: 



„2 



^4i(eJl7.|ei)(/R|7l/L) ,/R = iR ^^^^^ 

for top and bottom quarks. In all cases, log ^ terms can be savcly neglected to the 
accuracy we are working. Thus we find for left handed quarks of the third generation: 

A^Sro^p(Pi,-,P4) = MZM,-,Pd {l - lf^iS^W.log]i^} (140) 

For right handed external top quarks wc have 0^, x and H corrections. In that case 
we observe that the ZxH, 70=*" and Z4>^ loops have an opposite sign relative to the 
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left handed case. For the corrections corresponding to the topology shown on the right 
in Fig. 13 we must replace Qt in Eq. (138) by — T| = | for the 0=*" graph. The 
same contribution is obtained by adding the H and x loops and we find: 



At higher orders we note that the exchange of gauge bosons inside the one loop process 
is subsubleading and we arrive at the factorized form analogous to the Yukawa cor- 
rections for external Goldstone bosons. Since these corrections are of universal nature 
we can drop the specific reference to the Drell-Yan process and the application of the 
generalized Gribov-theorem for external fermion lines to all orders yields: 



where l^f^(s, M^) is given in Eq. (136) and the quantum numbers are those of the 
external fermion lines. Since at high energies all fermions can be considered massless 
we can again absorb the chiral top- Yukawa corrections into universal splitting functions 
as in Ref. [37]. Thus in the electroweak theory we find to next to leading order the 
corresponding probability for the emission of gauge bosons from chiral fermions subject 
to the cutoff fjL are given by Eqs. (128) and (130). The corrections from below the 
scale M need to be included via matching as described above. 

For physical observables, soft real photon emission must be taken into account in 
an inclusive (or semi inclusive) way and the parameter /^^ (134) will be replaced by 
parameters depending on the experimental requirements. This will be briefly discussed 
in section 3.4. Next we turn to longitudinal degrees of freedom after flrst reviewing 
the Goldstone boson equivalence theorem. 

3.2 The equivalence theorem 

At high energies, the longitudinal polarization states can be described with the 
polarization vector is given in Eq. (117). 

The connection between S-matrix elements and Goldstone bosons is provided by 
the equivalence theorem [105, 106, 107]. It states that at tree level for S-matrix ele- 
ments for longitudinal bosons at the high energy limit M'^/s — > can be expressed 
through matrix elements involving their associated would-be Goldstone bosons. We 
write schematically in case of a single gauge boson: 





M{pi, ...,p„;/x^) = Alnoop(Pi, ■■■,Pn) exp 




(142) 




(144) 



(143) 
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The problem with this statement of the equivalence theorem is that it holds only at 
tree level [108, 109]. For calculations at higher orders, additional terms enter which 
change Eqs. (143) and (144). 

Because of the gauge invariance of the physical theory and the associated BRST 
invariance, a modified version of Eqs. (143) and (144) can be derived [108] which reads 

^'^M(iy±(A;), Vphys) = C^M^At((/)±(^), V'phys) + O (145) 
A;'^M(Z,(A;), V'phys) = ^C.M,M (x(A;), V'phys) + O (146) 

where the multiplicative factors C„ and depend only on wave function renormal- 
ization constants and mass countcrterms. Thus, using the form of the longitudinal 
polarization vector of Eq. (117) we can write 

M{Wt{k),%^,,) = C^M{(j)^{k),ijp^,,) + O (^] (147) 



A^(Z4fc),Vys) =iC,M(x(fc),^phys) + 0(^^) (148) 

We see that in principle, there are logarithmic loop corrections to the tree level equiv- 
alence theorem. The important point in our approach, however, is that the correction 
coefficients are not functions of the energy variable s: 

= Cw(7Z, M, g, g') , C, = C,(7Z, M, g, g') (149) 

The pictorial form of the Goldstone boson equivalence theorem is depicted in Fig. 
14 for longitudinal W^-boson production at a linear e^e~ collider. In the following we 
denote the logarithmic variable t = log where is a cutoff on the transverse part 
of the exchanged virtual momenta k of all involved particles, i.e. 

Ii^<kl = mm{2{kpi){kpj)/{pipj)) (150) 

for all j ^ L The non-renormalization group part of the evolution equation at high 
energies is given on the invariant matrix element level by Eq. (84): 

^A^(L(A;),V^phys) = K{t)M{L{k),i^p^y,) (151) 

and thus, after inserting Eqs. (147), (148) we find that the same evolution equation 
also holds for M.{4>{k),^jJphys)■ The notation here is L = {W^, Z^} and = {0^,x}) 
respectively. Thus, the log ^ dependence in our approach is unrelated to the correc- 
tions to the equivalence theorem, and in general, is unrelated to two point functions 
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Figure 14: The pictorial Goldstone boson equivalence theorem for VF-pair production 
in e^er collisions. The correct DL-asymptotics for longitudinally polarized bosons arc 
obtained by using the quantum numbers of the charged would-be Goldstone scalars at 
high energies. 

in a covariant gauge at high energies where masses can be neglected. This is a conse- 
quence of the physical on-shell renormalization scheme where the MS renormalization 
scale parameter /I ~ M. Physically, this result can be understood by interpreting 
the correction terms C„ and as corrections required by the gauge invariance of the 
theory in order to obtain the correct renormalization group asymptotics of the physi- 
cal Standard Model fields. Thus, their origin is not related to Sudakov corrections. In 
other words, the results from the discussion of scalar QCD in section 2 should be apph- 
cable to the subleading scalar sector in the electroweak theory regarding a non-Abehan 
scalar gauge theory as the effective description in this range according to >Csymm in Eq. 
(11). The only additional complication is the presence of subleading Yukawa enhanced 
logarithmic corrections which will be discussed below. It is also worth noticing, that 
at one loop, the authors of Ref. [29] obtain the same result for the contributions from 
the terms of Eq. (149). In their approach, where all mass-singular terms are identified 
and the renormalization scale Ji = y^, these terms are canceled by additional correc- 
tions from mass and wave function counterterms. At higher orders it is then clear that 
corrections from two point functions are subsubleading in a covariant gauge. 

3.3 Longitudinal degrees of freedom 

According to the discussion of the previous section we can use Goldstone bosons 
in the high energy regime as the relevant degrees of freedom for longitudinal gauge 
boson production. Thus, the Higgs boson and the would-be Goldstone bosons actually 
receive the same corrections in high energy processes (up to purely electromagnetic 
terms). 
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Regulating the virtual infrared divergences with the transverse momentum cutoff 
as described above, we find the virtual contributions to the splitting functions for 
external Goldstone and Higgs bosons: 



Ti(Ti + l) + tan2^w (- 



-21og — 



+ 4 



3 

2M2 



5{1 - z) (152) 



The functions can be calculated directly from loop corrections to the elementary pro- 
cesses in analogy to QCD [65, 66, 67] and the logarithmic term corresponds to the 
leading kernel of Ref. [37]. We introduce virtual distribution functions which include 
only the effects of loop computations. These fulfill the Altarelli-Parisi equations in 
analogy to Eq. (74): 



d(t>{z,t) r^dy 



dt 



i: 



{z/y,t)PUy) 



(153) 



The sphtting functions are related by P^^ = + P^^, where R denotes the contri- 
bution from real boson emission. P^^ is free of logarithmic corrections and positive 

definite. 

Inserting the virtual probability of Eq. (152) into the Eq. (153) we find: 



4>{l,t) = 0oexp 



9' 



87r2 



Ti{Ti + 1) + tan^ 9^ 



3 m? s 



log' 



41og- 



2M2 



(154) 



These functions describe the total contribution for the emission of virtual particles 
(i.e. z = 1), with all invariants large compared to the cutoff /i, to the densities (f){z,t) 
{(f) = {0^, Xi H})- The normalization is not per line but on the level of the cross section. 
For the invariant matrix element involving external scalar particles we thus find at 
the subleading level: 



1 



M{pi,...,Pn,g,fi^) = AlBom(Pl,--,Pn,C/)exp|--^VI//(s,/i^ 



(155) 



where 



167r2 



Ti(Ti + l)+tanX 



^ V__41og-+-^log- 



(156) 



The functions Wf correspond to the probability of emitting a virtual soft and/or 
coUinear gauge boson from the particle 4> subject to the infrared cutoff ii. Typical 
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Figure 15: Two loop corrections involving Yukawa couplings of scalars to fermions. The 
Ward identity in Eq. (135) assures that in the Feynman gauge, the sum of both dia- 
grams does not lead to additional SL logarithms at the two loop level. Only corrections 
to the original one loop vertex need to be considered and lead to the exponentiation 
of Yukawa terms in the SM to SL accuracy. 



diagrams contributing to Eq. (156) in a covariant gauge are depicted in Fig. 8. The 
universality of the splitting functions is crucial in obtaining the above result. 

Again, since the Yukawa enhanced terms are novel features in broken gauge theories 
as compared to the situation in QCD we use the non-Abelian generahzation of the 
Gribov theorem in the following to prove the correctness of our splitting function 
approach for specific processes using on-shell renormalization of the external Goldstone 
bosons. 

Since the three fermion loop is more complicated than the situation in the fermionic 
sector above, we provide some more details in deriving the respective Ward identity. 
At the two loop level, we need to consider the diagrams displayed in Fig. 15. The 
corresponding two loop amplitudes read (neglecting I outside the fermion loop): 

J (477)^7 (47r)" {P - X^){k2 + iy(ki - lfk\k - kiYik -ki + If ^ ' 

r dn r d^k (fci - A;2).Tr [{GrUr + Gm){^ -h + ^(^^r + Giuji)^ Ak^k^ 
J {An)^] (47r)" (P - X^){k2 + lf{ki - Ifk^ik - ki + lY{ki - If ^ ' 

where we omit common factors and the scalar masses taking M ~ A for clarity. The 
soft photon corrections must also be included via matching. The G^,/ denote the chiral 
Yukawa couplings and u}r,i = | (1 ± 75)- The gauge coupling is again written in the 
symmetric basis. For our purposes we need to investigate terms containing three large 
logarithms in those diagrams. Since the fermion loops at one loop only yield a single 
logarithm it is again clear that the gauge boson loop momentum / must be soft. Thus 
we need to show that the UV logarithm originating from the k integration is identical 
(up to the sign) in both diagrams. We can therefore neglect the loop momentum / 
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inside the fermion loop. We find for the fermion loop vertex r^{kl, 0, kj) belonging to 
Eq. (157): 

_ AGrGi {2kt{k'^ - kik) + k^jkl - k^)) 
~ k^{k - kiY 

This we need to compare with the self energy loop from Eq. (158): 

d Tr [{GrUr + GiUi){1^ - ^^){GrUr + GiUi)^ 



(159) 



2 



9A;i^ k'^{k — ki) 

d AGrGi{k^k-e) 2k^,{k^-hk) + k^{kl-k^) 



dkx^ k^ik-ki)"^ " ' k^ik-kiY 

In short we have estabhshed the analogous Ward identity in Eq. (135) in the lon- 
gitudinal would-bc Goldstone boson sector. The relative sign is again such that the 
generated SL logarithms of the diagrams in Fig. 15 cancel each other. The existence of 
such an identity expresses the fact that also the longitudinal Yukawa sector is gauge in- 
variant. At higher orders this Ward identity ensures the corresponding cancellations to 
SL accuracy. Also in an axial gauge the corrections can be seen to factorize accordingly 
since in this gauge DL terms originate only from on-shell two point functions. 

We are thus left with gauge boson corrections to the original vertices in the on-shell 
renormalization scheme such as depicted in Fig. 16. At high energies we can therefore 
employ the non-Abelian version of Gribov's bremsstrahlung theorem. 

Thus the corrections at higher orders factorize again with respect to the one loop 
amplitude and M. "Bom" {Pi, ■■■,Pn) = -Miioopipi, ■■■■iPn)- Note that the latter is also 
independent of the cutoff fi since the fermion mass serves as a natural regulator. The 
subleading corrections are then compared to the general splitting function approach 
yielding Eq. (155). 

In the case of the amplitude of Fig. 16 we must use the quantum numbers of the 
associated Goldstone bosons and we have the following Born amplitude 

A^Bom(pi, ■■■,P4) = ^^^(eR|7''|eL)(^i - ^2)1. (161) 

and at one loop we have two fermion loops contributing {ttb and bbt). The renormaliza- 
tion condition is provided by the requirement that the corrections vanish at the weak 
scale, i.e. for s = M^, which amounts to subtracting the vertex for that case. The 
first diagram of the two is, according to the Feynman rules of appendix 6.2 using the 
physical fields, given by 



4 2 



+ x 
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Figure 16: A Feynman diagram yielding Yukawa enhanced logarithmic corrections in 
the on-shell scheme. At higher orders, the subleading corrections are given in factorized 
form according to the non-Abelian generalization of Gribov's theorem as described in 
the text. Corrections from gauge bosons inside the top-loop give only sub-subleading 
contributions. DL-corrections at two and higher loop order are given by gauge bosons 
coupling to (in principle all) external legs as schematically indicated. 



/ 



(27r)" (Z2 - ml + %e){{l + k^Y - m\ + %e){{l - k^Y - m\ + is) 

3iQt e'^ml „m 



{e^\Y\et){B,,-B^,){k,-k,), 



(162) 



where (jJ± = |(1±75) and the chiral couplings are given by c£ = Qf for the photon and 

c+ = ^Qf ^-iid c_ = ^ ^ ^ — ^ for Z-bosons respectively. The counterterm 5*f is chosen 
such that the logarithmic corrections vanish for s — M^. Thus, the sum of the scalar 
functions is to logarithmic accuracy S23 — — —\og^. Analogously, we have for 
the hht quark loop: 



/ 



Hoop 



3E 



4 2 
e ml 

2sM^sl 



hht 



{2tiY (/2 - ml + ie){{l + kiY - ml + ie){{l 



k2Y -ml + is) 



+ « 



leTT^ci 



4 2 

e mi 



{e^\Y\el){B2^-B^^){k,-k2). 



(163) 
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H , 

Figure 17: A Feynman diagram yielding Yukawa enhanced logarithmic corrections to 
external longitudinal Z-bosons and Higgs lines in the on-shell scheme. At higher orders, 
the subleading corrections are given in factorized form according to the non-Abelian 
generalization of Gribov's theorem as described in the text. Corrections from gauge 
bosons inside the top-loop give only sub-subleading contributions. 



Adding both results (162) and (163) we find 

{g^ 3 Tfl^ S I 

^ ~ T6^2M^^°^M2 J ^^^^^ 
and the all orders result to subleading accuracy is given by 

M(pi, A^') = Mnoop(Pi, -,Pn) exp (^-^ Wr{s, M^) j (165) 

where Wf'^{s^ M^) is given in Eq. (136). The subleading Yukawa corrections from the 
Altarelli-Parisi in Eq. (155) agree with the corresponding results from the application 
of the Gribov-theorem in Eq. (165). For longitudinal Z-boson and Higgs production, 
we note that there is only one non-mass suppressed elementary vertex with two neutral 
scalars, namely the ZxH vertex. As mentioned above, universal terms are related to 
the massless limit. For the "Born amplitude" of the Higgs-strahlung vertex we have 

The universal Yukawa corrections to both external x ^-^d H states from an off shell Z 
line are then given by the corrections depicted in the inner fermion loop of Fig. 17. 
Here we find 



^f4>,-,P3)=3^x 
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(27r)" (Z2 - m? + ie){{l + k^Y - mf + ie){{l - - mf + ie) ^ 

fi7^3 3 2 

* *-(i?23-i?^)(fcl-fc2) (167) 



and thus 



M^CiP^^ ...,Ps) = M^^f. ] 1 - Iy^^, log ^ \ (168) 



Prom the same hne of reasoning as for the charged Goldstone bosons we find that the 
all orders result is given by Eq. (165). At the subleading level, this is equivalent to 
the corresponding corrections obtained in Eq. (155). 

In addition to the Sudakov corrections in Eq. (156) we also have to include terms 
corresponding to the renormalization of the mass terms in the Yukawa coupling of 

the Born amplitude 

^ 2 5 j ^t the one loop level [29]. At higher orders, mass 

renormalization terms are connected to two point functions and thus subsubleading. 
The relevant higher order SL-RG terms, however, will be discussed section 3.5. 

3.4 Semi-inclusive cross sections 

In order to make predictions for observable cross sections, the unphysical infrared 
cutoff /i^ has to be replaced with a cutoff /^expf related to the lower bound of fc^ 
of the other virtual particles of those gauge bosons emitted in the process which are 
not included in the cross section. We assume that A*expt < that the non- 

Abelian component of the photon is not essential. The case /igxpt > much more 

complicated and is discussed in Ref. [28] through two loops at the DL level. 

We can write the expression for the semi-inclusive cross section in a compact way 
as follows: 

da(pi, ...,pn,g, g', //expt) = daeiasudPi, ...,Pn,9, g', 

X exp('u;Jxpt(s, mi, ji, /^expt)) (169) 

The dependence in this expression cancels and the semi-inclusive cross section de- 
pends only on the parameters of the experimental requirements. In terms of all the 
virtual corrections discussed above the full result is given by: 

da{pi, ...,pn,g, g', /^expt) = daBomiPi, ...,Pn,g, g') 

( "g "/ ^ 

X exp - ^ W,, {s, M')-Y: Wf, {s, M')-Y: W^, {s, M') 



X exp 



i=l i=l i=l 



1=1 



1=1 



X exp (w^expt(s> rui, //expt)) (170) 
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where Ug denotes the number of transversely polarized gauge bosons and rtf the number 
of external fermions. This expression omits all RG corrections, even at the one loop 
level. The functions W and w correspond to the logarithmic probability to emit a soft 
and/or collinear particle per line, where the capital letters denote the probability in 
the high energy effective theory and the lower case letter the corresponding one from 
pure QED corrections below the weak scale. The matching condition is implemented 
such that for n = M only the high energy solution remains. For the contribution from 
scalar fields = {0^, Xt H} above the scale M we have W^^{s^ M^) in Eq. (156) with 
a = g'^/4:TT and tan^ 6*^ = a' /a. The last term is written as a logarithm containing 
the top quark mass rrit rather than the weak scale M since these terms always contain 
nit as the heaviest mass in the loop correction [38]. M^) and Wg. are given in 

Eqs. (128) and (129) respectively. Again we note that for external photon and Z-boson 
states we must include the mixing appropriately as discussed in section 3.1. For the 
terms entering from contributions below the weak scale we have W/.(s, /x^) given in 
Eq. (131), in Eq. (132) and Wy^ in Eq. (133) for the virtual corrections. For real 
photon emission we have in the soft photon approximation: 



^7expt(s>^i>/^>/^expt)= < 



^i=l (47r)2 
l^i=l (47r)2 



-log^ 



-'cxpt 



+ log2-^-31og. 



log A - 1 
-2 log 



21og^ + log2 



,mi <^ fj, 



(171) 



log 



where n is the number of external lines and the upper case applies only to fermions since 
for we have fx < M. Note that in all contributions from the regime jj, < M we have 
kept mass terms inside the logarithms. This approach is valid in the entire Standard 
Model up to terms of order O (log . Note, however, that in the Yukawa enhanced 
terms the replacement M — > rrit is legitimate. The overall yU-dependence in the 
semi-inclusive cross section cancels and we only have a dependence on the parameter 
A^expt related to the experimental energy resolution. All universal electroweak Sudakov 
corrections at DL and SL level exponentiate. 

3.5 Electroweak RG corrections 



The way to implement the SL-RG corrections is clear from the discussion in section 
2.6. At high energies, the DL phase space is essentially described by an unbroken 
SU{2) X U{1) theory in which we can calculate the high energy contributions. In this 
regime, all particle masses can be neglected and we have to consider the following 
virtual electroweak DL phase space integral with running couplings in each gauge 
group: 
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27r [c \ / q;(s) J 

where Q;(/i^) = g'^ {/j,'^) / An , a!{p?) — g''^(//^)/47r, c = a{iJ,'^)f3o/7r and analogously, c' = 
a' {fi'^)j3'Q/7i. In each case, the correct non-Abelian or Abelian limit is reproduced by 
letting the corresponding couplings of the other gauge group approach zero. In this 
way it is easy to see that the argument of the running couplings can only be what 
appears in Eq. (172). 

The form of Eq. (172) is valid for fermions, transversely and longitudinally po- 
larized external lines but (omitted) subleading terms as well as the quantum numbers 
of the weak isospin Tj and the weak hypercharge Yi differ. In order to implement 
the missing soft photon contribution, we choose the analogous form of solution in Eq. 
(114) and have to implement it in such a way that ior /j, — M Eq. (172) is obtained. 
The full result for the respective semi-inclusive cross sections is then given by: 

dcr^^ipi, ■ ■ ■ ,Pn,9,9',fJ'cxpt) = daBorniPi, ■ ■ . ,Pn,g{s),g\s)) 

( "9 ^ 

X exp - ^ W^^is, M') - ^ Wf^{s, M') - Wf^{s, M') 
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-^^,,(M^m^2) 



1=1 



exp (w^^^^ {s, rrii, /Xexpt)) (173) 



where Uf denotes here again the number of external fermions. The argument of the 
gauge couplings in the Born cross section indicate the one loop renormalization of the 
couplings which is not included in the exponential expressions but which at one loop 
is genuinely subleading: 

a{s) = a{M') (l - Po^^ log (174) 



TT 



a' is) = a'{M^) (l - P',^^ log ^) (175) 
where a{M^) = e'^{M^)/A7rsl and ^'(M^) = e'^{M'^)/A7rcl with 

Am') = (i + ~ E Qn i«« (176) 

and e'^/AiT = 1/137. If there are non-suppressed mass ratios in the Born term, also 
these terms need to be renormalized at one loop (see Ref. [29]). Higher order mass 
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renormalization terms would then be sub-subleading. The function W^^{s, M^) is 
given by 
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where we again have rrit in the argument of the Yukawa enhanced correction 
Analogously for fermions we have: 
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The last term contributes only for left handed bottom and for top quarks as mentioned 
above and /' denotes the corresponding isospin partner for left handed fermions. 
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Again we note that for external photon and Z-boson states we must include the mixing 
appropriately as discussed in Ref . [37] . For the terms entering from contributions below 
the weak scale we have for fermions: 
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where c = — <5j-/V^- Analogously, for external W-bosons and photons we 

find: 
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At < 

Note that the function w^^{M'^,ii^) does not receive any RG corrections to the order 
we are working since it contains only SL terms. For the virtual corrections and for real 
photon emission we have in the soft photon approximation: 
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where n is the number of external lines and Uf fermions propagating in the loops folded 
with the DL integrals. The upper case applies only to fermions since for we have 
ji < M. Note that in all contributions from the regime /j, < M we have kept mass 
terms inside the logarithms. For the running above the weak scale M we use only the 
massless Po, P'q terms with Ugen = 3. This approach is valid in the entire Standard 
Model up to terms of order O (\og^^. 

In order to briefly discuss the size of the SL-RG corrections obtained here, we 
present a numerical comparison. For this purpose we will only compare the terms 
which are new in the RG analysis, i.e. the running from the weak scale M to y/s. 
We are thus interested in effects starting at the two loop level and want to compare 
the relative size of the RG-improved form factors to the pure Sudakov terms. It is 
therefore of interest to compare the ratios (e^~^i^'^^ — e^^^^-^j /e^~^i^^^ for the various 



particle labels i. Since the physical scales in the problem are given by M and ^/s, the 
lower and upper limits of the couplings are given accordingly by these scales for the 
functions Wj. Fig. 18 compares the respective ratios for various SM particles. For 
definiteness we take M = 80 GeV, rut = 174 GeV, si = 0.23, a{M^) = 1/128/s^, 
a'{M^) = 1/128/4, Po = 19/24 and = -41/24. The difference between the curves 
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using and those using s as the scales in the conventional Sudakov form factors 
is a measure of the inherent scale uncertainty which is removed by the RG-improved 
Sudakov form factors W^^. The largest effect is obtained in the gauge boson sector 
[77]. For external j^"*", %, i7} particles we have at 1 TeV a difference between the 
curves of about 0.35% per line on the level of the cross section, growing to 0.65% 
at 2 TeV. The situation is very similar for transversely polarized W'^, W' particles 
where it reaches about 0.4% at 1 TeV and 0.8% at 2 TeV per line on the cross section 
level. For left handed quarks of the third generation the size of the corrections is about 
0.15% at 1 TeV per line on the level of the cross section and 0.33% at 2 TeV. These 
corrections are thus considerably smaller and only needed if precisions below the one 
percent level are necessary from the theory side. For right handed top quarks the effect 
is even smaller since only the running of a' enters and it is thus negligible for most 
applications. The form of the two curves in case of right handed tops differs markedly 
from the other three cases because at the energies displayed, the dominant effect is 
actually due to subleading Yukawa enhanced corrections (~ a) since the DL terms are 
proportional to a' and since the ratio ml/M'^ is of the size of an additional logarithm 
for these values of 

In general it can be seen that-where the DL terms dominate-the renormalization 
group improved results are indeed in-between the upper and lower bounds given by 
the respective scale choices in the conventional Sudakov form factors. Indeed also for 
right handed top quarks this pattern is observed if only DL corrections are taken into 
account. 

It should be emphasized again that also the QED-RG corrections can be sizable 
since large mass ratios with light particles occur. These should of course also be 
implemented in a full SM prediction at TeV energies keeping in mind that one must 
always be far above any particle thresholds. 

3.6 Angular dependent corrections 

In this section we discuss the important contribution from angular dependent log- 
arithmic corrections. In massless four fermion processes these were first given by Ref. 
[27] at the SL level and in Ref. [110] at the SSL level. Since the size of the observed 
corrections is large, these terms are important for future collider phenomenology. In 
Ref. [Ill] a general way of treating SL angular terms for arbitrary processes including 
mass terms was presented. The regime of validity assumes that all 2piPj S> M^. In or- 
der to better understand the origin of angular dependent corrections in the electroweak 
theory, consider the case of massive QED and the right handed SM depicted in Fig. 
19. We will show that the mass terms do not lead to new effects compared to those in 
the massless case. 
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(expi-wf'^} - exp{-Wi}) / exp{-wp°} 
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Figure 18: This figure compares the renormahzation group improved probabihties Wj 
with the conventional Sudakov exponentials Wi for various external particle lines. The 
comparison is made with the indicated scale choices for the functions Wi and takes 
into account only the RG corrections from the scale M to Taking the difference 
between the two curves is a measure of the uncertainty removed by the RG effects. 
The variations in the scale of the couphng in the Wi functions is largest in the scalar 
(Goldstone and Higgs boson) sector and for transverse W^ where the effect is about 
0.8% at 2 TeV per line on the level of the cross section. In general, the RG improved 
form factors differ by fractions of one percent per line and need to be taken into account 
at future colliders if the experimental accuracy is in the percentile regime. 
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P3 Pa 




Pa Pa 



Figure 19: Angular dependent two loop on-shell diagrams. The sum of 11^ and 11^ 
factorizes in massive QED and the right handed SM into the product of the two one 
loop corrections (each with a different invariant and mass terms) in leading order. 



The two scalar integrals of Fig. 19, regularized with gauge boson mass terms Ai 
and A2, are given in massive QED or in the case of right handed (massive) fermions by 

ir = Astl^(^ { ^ X 

' J (27r)4 J (27r)4 \ (Ij - Af)((pi - h)^ - mj)((p, -h- hf - ml) 

(185) 



{{P2 + hr - ml){ll - Ai)((p3 + kf - ml) 
' J (27r)4 J (27r)4 \ {If - A?)((pi - h)' - m?)((pi -h- 

{{p, + h)^ - ml){ll - Ai)((p3 + k)' - ml) } ^^^^^ 

denoting s = 2piP2, t — 2piP3 and where the m^ are the masses of the external charged 
particles on their mass shell. Thus, it is straightforward to see that the sum of the two 
diagrams factorizes to leading order: 

ir + tf = 4.t/^ f^i il + il-2p,{h + k) ^ 

' ''J {2n)^J {2n)^\{q-Xl){{p,-l,)2-ml)((p,-h-k)'-m'i) 
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The omitted cross term 2/1/2 leads only to corrections containing three logarithms at 
the two loop level. It is thus on the same level as the approximation in the beginning of 
our discussion which only considers scalar integrals and can therefore be neglected. To 
DL accuracy we can employ the Sudakov technique, parametrizing the loop momenta 
along the external four momenta as 
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Thus, after rewriting the measure and integrating over the perpendicular components 
we find for the case of two photons with mass A (omitting the principle value parts): 
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The important point about the result in Eq. (190) is not only the factorized form in 
terms of the two massive one loop form factors but also the fact that the fcrmion mass 
terms correspond to each external on shell line in the amplitude. Thus, by rewriting 
the term in the bracket of the last line in Eq. (190) as 
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we see that the SL angular terms are indeed independent of the fermion mass terms. 
For the corrections involving the invariant u = 2p2P3 the situation is analogous. 

In the case of the SM with right handed massive fermions we need to consider in 
addition the exchange of Z-bosons. The results read 
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for the case of a photon and a Z-boson with mass M, and for the case of two Z's we 
have 
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Again we see the independence of the SL angular terms on the fermion mass terms 
and in addition, the fact that the gauge boson mass gap does not spoil the type of 
factorization in the right handed SM. 

This type of factorization can be generalized on theoretical grounds to the situation 
in the general SM. It should be noted that all fermion mass singularities in the SM only 
arise through photon radiation or coupling renormalization. The latter is not important 
in our discussion here and is anyhow sub-subleading at higher orders. The exchange 
of the heavy gauge bosons does not lead to fermion mass singular terms assuming that 
all rrii < M. This case is analogous to QCD where angular terms factorize in matrix 
form [86]. Thus, only corrections where one heavy gauge boson and one photon are 
involved are novel features in the SM. In this case, however, the type of factorization 
analogous to Eq. (190) for fermion mass and soft terms follows from the factorized 
form of real emission corrections. They are of factorized form and the KLN theorem 
then leads to the analogous situation for the sum of all virtual corrections. The soft 
terms must also factorize since we can always define an observable by only allowing 
soft photon radiation, i.e. chosing < M. 
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Prom the arguments presented in Ref. [Ill] it then follows that 
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where denotes the number of external gauge bosons (in the symmetric basis), n/ the 
number of external fermions and the number of external scalars (including Higgs 
particles). The notation used in the regime below the scale M is analogous; note that 
we use a photon mass regulator in this case. The fields u have a well defined isospin, 
but for angular dependent terms involving CKM mixing effects, one has to include 
the extended isospin mixing appropriately in the corresponding couplings TT". of the 

k ' ^ 

symmetric basis. The expressions for the virtual probabilities W^'^ are given in section 
2.6. In addition we denote 



for the angular dependent corrections from the soft QED regime, where the Qj denote 
the electrical charges of the external lines. These terms are not manifestly positive and 
thus not related to emission probabilities and also process dependent. 

The diagonal terms do not involve (CKM-extended) isospin rotated Born matrix 
elements. These occur only from the angular terms above the scale M. Thus, the way 
the electroweak angular terms factorize is also in an exponentiated operator form as 
in QCD. Eq. (194) is valid for arbitrary external lines. 

3.7 Sudakov logarithms in softly broken supersymmetric models 

In the introduction we already discussed the motivations for considering supersym- 
metric extensions of the SM. Since supersymmetry must be broken in nature (if it is 
at all relevant to physics at the electroweak scale) , the original simplicity is lost and in 
the case of the MSSM over one-hundred parameters are needed to describe the model. 
In order for supersymmetry to stabilize the hierarchy problem, though, two conditions 
need to be fulfilled. Firstly, the masses of the superpartner particles must not be much 
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larger than the electroweak scale and secondly, it must be a "softly broken" symmetry, 
i.e. broken by mass terms and couplings with a positive mass dimension. 

Thus, at energies in the TeV regime, these two conditions (assumed to be fulfilled in 
the MSSM or the NMSSM for instance) lead to the following consequences for radiative 
corrections of electroweak origin. If we assume that the mass scale of the superpartner 
particles are not much larger than the weak scale, say less than 500 GeV, and energies 
in the TeV regime, the results for the DL corrections outlined above for the SM can 
be applied straightforwardly to the MSSM. The reason is that the gauge couplings are 
preserved under supersymmetry and no additional spin 1 particles are exchanged. The 
appropriate quantum numbers in the eigenvalues of the casimir operators are the same 
as those of the SM partners. 

In case the superpartner masses are larger than 500 GeV, additional double loga- 
rithms need to be taken into account in a way outlined in Ref. [28]. In the following 
we assume that we can neglect such terms, i.e. that all particles in the MSSM have a 
mass below 500 GeV. 

At the subleading level, the situation in general is less clear at higher orders. For 
SL angular dependent terms, the same reasoning as above goes through since they 
originate only from the exchange of spin 1 gauge bosons and can thus be resummed as 
in the SM (see section 3.6). Box-type diagrams exchanging supersymmetric particles in 
the s-channel do not contribute to SL angular terms. The same holds for all universal 
SL corrections which involve the exchange of SM particles since they are properties of 
the external particles only. 

New types of SL Sudakov corrections are, however, involved in the exchange of 
supersymmetric particles as discussed in Refs. [112, 113, 114] at the one loop level in 
the on-shell production of superpartner particles in e'^e" collisions. In the following we 
discuss the results obtained in those works. We begin with the corrections contributing 
in particular the Yukawa terms from the final state corrections. The final result of 
calculating all terms contributing to the Yukawa sector includes terms depending on 
tan/3, which is the ratio of the two v.e.v.'s from the two Higgs doublet sector of the 
MSSM. Since tan/3 could be as large as 40, the bottom Yukawa terms are crucial in 
supersymmetric models. Also higher order terms need to be considered as in the SM 
and the exponentiation of the MSSM Yukawa terms follows from the same arguments 
as in the SM. In particular the Ward identity Eq. (135) holds [114]. Again, here 
we assume that the susy masses are close to the electroweak scale. These corrections 
apply to both the SM as well as the superpartner production. In particular the process 
of charged Higgs production seems to be well suited for an indirect measurement of 
tan^ > 14 of better than (9(25%), and a few percent for tan/3 > 25 [114]. The 
important point to note here is not only the precision but in particular the fact that 
this determination of tan /3 is independent of soft breaking terms to SL accuracy and 
thus model independent. This is due to the fact that the soft breaking contributions 
are constants and can be eliminated via subtraction if a series of precise measurements 
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is performed at various energy scales. In addition this approach is scheme and gauge 
invariant. 

Also the SL-gauge terms get modified by loops containing novel superpartner con- 
tributions. If we neglect the mass difference between the SM fermionic particle and its 
scalar partner we observe an exact supersymmetry relation^ described by (for A = M) : 
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where Tj denotes the total weak isospin of the particle j, Yj its weak hypercharge and 
at high s the invariants are given by t = — | (1 — cos^) and u = — | (1 + cos 6*). The 
helicities are those of the fermions (/) whose superpartner is produced. In addition 
we denote rhj — rrit/ sin/3 if / = f and rhj — m^j cos/3 \i j — h. f denotes the 
corresponding isopartner of /. For particles other than those belonging to the third 
family of quarks/squarks, the Yukawa terms are neghgible. Eq. (196) depends on the 

^Here we take a common mass also for the two chiral superpartners. Since in general the mass 
eigcnstatcs arc different from the flavor cigcnstatcs if their masses differ when supersymmetry breaking 
is included, mixing effects must then also be considered. 
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important parameter tan (3 — the ratio of the two vacuum expectation values, and 

displays an exact supersymmetry in the sense that the same corrections are obtained 
for the fermionic sector in the regime above the electroweak scale M. 

Here we assume that the asymptotic MSSM /5-functions can be used from the scale 
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where Ca = 2, Ug = 3 and = 2. In practice, one has to use the relevant numbers 
of active particles in the loops. These terms correspond to the RG-SL corrections just 
as in the case of the SM as discussed in Ref. [77] but now with the MSSM particle 
spectrum contributing. They originate only from RG terms within loops which without 
the RG contribution would give a DL correction. It should be noted that the one-loop 
RG corrections do not exponentiate and are omitted in the above expressions. They 
are, however, completely determined by the renormalization group in softly broken 
super symmetric theories such as the MSSM and sub-subleading at higher than one 
loop order. 

The relation (196) is expected since in unbroken supersymmetry both chiral 
fermions and the superpartner sfermions are part of the same supermultiplet. In the di- 
agrammatic evaluation, however, the identity expressed in (196) is a highly non-trivial 
check on the overall correctness of the calculation as different particles and loops con- 
tribute in each case. Since in the real world ^ ruf, the corresponding matching 
terms containing light fermion mass need to be included as well. 

For charged Higgs production we have analogously 
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It should be noted here that the Yukawa terms proportional to tan/? are quite large 
due to the additional factor of 3 = Nc from the quark loops [114]. Overall, both 
Yukawa contributions in Eq. (196) and (199) reinforce the Sudakov suppression factor 
of the leading double logarithmic terms. 

Eq. (199) is also valid for charged Higgsino production (as suggested by the relation 
(196)), however, the latter is not a mass eigenstate. 

The higher order exponentiation of the universal terms uses the Ward identity 
(135) for the respective corrections. The angular terms are of SM origin only as de- 
scribed above and are thus treated as in section 3.6. In addition we have the matching 
terms stemming form mass terms of both the gauge bosons and particles involved in 
the process. In any case it would also be helpful for collider experiments to have a 
full one loop calculation in the full supcrsymmetric theory in order to have a better 
understanding of the size of constant terms omitted in our high energy approximation. 
These terms, however, are difficult to discuss at this point since they depend strongly 
on how supersymmetry breaking is realized in nature. 

Neutral Higgs production in the MSSM depends also on the parameter a, the 
mixing angle between the CP even neutral Higgs particles. For these terms the cor- 
responding dependence on tan/3 is thus not as evident as in the charged Higgs case 
without independent knowledge on a but results can be found in Ref. [114]. For trans- 
verse gauge boson production it should be clear that the (3 functions in Eqs. (179) 
should be replaced by the corresponding MSSM expression given in Eq. (197). It 
can also be expected that an analogous supersymmetry relation should hold between 
the high energy splitting function expression of gauge bosons and gauginos under the 
above assumptions (in particular the neglect of the mass differences) . 

Longitudinal gauge boson production is analogous to Higgs production after the 
application of the Goldstone equivalence theorem. 

3.8 Fully-inclusive cross sections 

Physical observables in the electroweak theory depend on the infrared cutoff, the 
gauge boson mass M, as was discussed in section 3.4. This feature is closely related 
to electroweak symmetry breaking itself since it indicates that in the unbroken limit 
M — > cross sections become unobservable. Thus, we would expect in this limit 
confining effects limiting the types of initial or final state particles to be in an isospin 
singlet state. In the broken physical theory, this is evidently not the case and physical 



71 



particles carry a non-Abelian group charge. It is for this reason that even fully inclusive 
cross sections are expected to depend on log^ Also in QCD the Bloch-Nordsieck 
violating terms are present in fully inclusive cross sections and only the initial state 
color averaging eliminates the infrared problem. Thus we see that masslessness of the 
gauge bosons in a non-Abelian theory is intimately connected with confinement and 
therefore color neutral initial and final states. 

The realization of this phenomenon in the electroweak theory was emphasized in 
Ref. [33] where these "Bloch-Nordsieck violating" terms were calculated by means of 
the coherent state operator formalism from QCD. The basic idea is provided by the 
observation that in collider experiments one does not take the average over, say, the 
e~e^ and ve^ initial states since the asymptotic states carry weak isospin quantum 
numbers. The authors of Ref. [115] conclude that these log^ terms in fully inclu- 
sive electroweak processes are due only from initial state corrections. There are, 
however, possibly additional terms from the Abelian sector due to the fact that the 
hypercharge is broken in the physical fields. These terms, however, can be neglected 
for light fermions. By then considering the leading coherent state operators for the 
soft and hard (Born) parts of the scattering amplitude, the authors of Ref. [115] find 
that the Bloch-Nordsieck violating terms exponentiate in the form of the Sudakov 
form factor in the adjoint representation. For initial state fermions with two isospin 
components the results for the fully inclusive cross sections are then given by 
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is the Sudakov form factor and the factor of 2 in the exponential argument corresponds 
to the Casimir eigenvalue in the adjoint representation, i.e. T(T + 1) with T = 1. A 
generalization to transversely polarized initial states is given in Ref. [116]. Denoting 
the triplet representation as -|-, 3, — the authors obtain for the case of fermion boson 
scattering: 
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In the case of transverse boson boson scattering we have 



(203) 
(204) 
(205) 



a 



Born 
'++ 



^exp 



:-2^.s 



1 

-exp 

D 



-6^5) + (T 



,Born 



1 _ 1 exp (-2^5) 



72 



+ 1 exp (-6^5)) + <jfr {I -I exp (-6^5)) (206) 
<7_ + = a+ _ = <7^r {I -I exp (-2^5) + I exp (-6^5)) + a^T + ^ exp (-2^5) 



+ - exp (-6^5) ) + ( o - o exp (-6^5) ) (207) 



.iexp(-6^.))+.3rQ4 

(73+ =a3_ = ^{[a^r + <^^T](l-exp(-6^5)) + ^T3^r(l + 2exp(-6^5))} (208) 

In Ref. [117] results for the longitudinal and Higgs sector are presented employing 
the equivalence theorem. For a light Higgs the inclusive cross sections are related as 
follows: 

(^(l>-<t>- = cr4>+4)+ ; cr4>+h = cr<p-h = cr<p+3 (209) 
CThh = (733 = ^ {cr^+4>+ + cr4,+^-) + Re (7^°^*^^.^^ exp (-(2 + tan^ 9^)J^s) (210) 

cr^h = ^ ((70+<^+ + t7,/>+0- ) - Re (tJo°^?^<^20o exp (-(2 + tan^ ^w)^^) (211) 

The individual results can be obtained from 

cr</,+</,+ + (^<p-4>+ + 2cr3^+ = o-^+^+ + 0-^°™+ + 2af^r (212) 

+ a^-4>+ - 2^3^+ = + aJ-7+ - 2af;r) exp (-2^^) (213) 

o-^+</,+ - cr^-4>+ + "^Im (T^4>+-^h<t>+ = <^^+^+ - + 2/m o-£^°T'^fe<^+ (214) 

- <7<^-^+ - 27m (J3^+_,<^+ = - <T?-7+ - 2/m <7f;r_,^+) exp (-2^5)(215) 

In general, the effects are sizable. In particular for initial state leptons and for gauge 
bosons O (10%) effects can be expected in the TeV range. In Ref. [118] an initial 
step towards electroweak splitting functions is taken for a spontaneously broken SU (2) 
gauge group. In this model there is no gauge boson mass gap, however, the important 
difference to the QCD case is that the initial state isospin averaging is omitted. A 
rigorous extension to the full electroweak theory would be important since the size of 
SL terms in non-negligible. 

3.9 Physical fields at fixed order 

In this section we present results obtained in real fixed order pcrturbative calcu- 
lations performed with the physical SM fields in the on-shell scheme. It is pivotal to 
compare these results with the corrections obtained with the IREE method expounded 
on above where the calculation is performed in the high energy effective theory includ- 
ing appropriate matching conditions at the weak scale in order to include pure QED 
effects. Fortunately, a lot of effort has been spent in pcrturbative electroweak calcula- 
tions at high energies in order to allow for a cross check for all degrees of freedom at 
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one loop. At the two loop level calculations at the DL level are available in the regime 

we are interested in. We begin, however, at the lowest order. It should be mentioned 
that in order to obtain an accuracy in the 1 % regime, it is not sufficient to include 
only logarithmic corrections since for processes involving unsuppressed mass-ratios in 
the Born cross section (such as — > 4>'^4>~) constants can be much larger. Also 
TT^ terms and Yukawa enhanced constants can be of the order of several % and thus, 
for the precision objectives of the linear collider a full one loop calculation is necessary 
in most cases. At higher orders, however, these constant terms can be neglected. 



3.9.1 One loop results 

There have been various calculations at the one loop level, however, for our pur- 
poses it is convenient to discuss the results of Ref. [29] in more detail where general 
formulas are derived for one loop logarithmic corrections at high energies. For all SL 
Sudakov logarithms we find agreement with the results presented in Ref. [29]. The 
universality of the corrections in the physical theory, however, is more difficult to see 
since the Ward identities are more involved in the broken gauge theory. We summarize 
results presented in Ref. [29] for DL and one loop angular terms containing log log | 
type contributions. In addition there are also the on-shell parameter renormalization 
of the couplings, masses and mixing angles. All results are in the 't Hooft-Feynman 
gauge and the techniques discussed in this section are valid only at the one loop level. 

In order to obtain double logarithmic corrections one can drop the dependence of 
the numerator on the loop momentum. In this approximation the one-loop corrections 
give 
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and in leading order, using the high-energy expansion of the scalar three-point function 
], one obtains 
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where the short hand notations 

i.(h,|.M^):^^W^ i(.«.M^):^^logl5L (218) 
L{s} ■-L(s,M^) 1(3} --Us, Mi,) (219) 

are used with 

Tki - {Pk + Pif ~ '2pkPi » (220) 
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The DL term containing the invariant r^i depends on the angle between the momenta 
Pk and pi- Writing 

L{\ni\,M^) = L{s, + 2/(s, M2) log ^ + L{\ruls) (221) 

the angular-dependent part is isolated in logarithms of rfc//s, and gives a subleading 
soft-coUinear (SSC) contribution of order l{s) log(|rfc;|/s), whereas terms L(|rfc/|, s) can 
be neglected in LA. The remaining part, together with the additional contributions 
from photon loops in (217), gives the leading soft-collincar (LSC) contribution and is 
angular-independent. This contribution corresponds to the universal part discussed in 
the previous sections in the symmetric basis. The eikonal approximation (216) applies 
to chiral fermions, Higgs bosons, and transverse gauge bosons, and depends on their 
gauge couplings /^"(/c) which can be found in appendix 6.1. 

Since the longitudinal polarization vectors (117) grow with energy, matrix elements 
involving longitudinal gauge bosons have to be treated with the equivalence theorem, 
i.e. they have to be expressed by matrix elements involving the corresponding Gold- 
stone bosons. As discussed in section 3.2, the equivalence theorem for Born matrix 
elements receives no DL one-loop corrections. Therefore, the soft-collinear corrections 
for external longitudinal gauge bosons can be obtained using the simple relations 



rDL A yl-'Z.T.... _ ^rDL 



^UL_^...Zl... ^ -^DL^^.^x... (222) 

from the corrections (217) for external Goldstone bosons. 

The invariance of the S matrix with respect to global SU(2) x U(l) transformations 
implies (up to mass terms): 

= ^v-^A^^i-^" = ie V lY^. {k)M''-''''-'- (223) 

k 

which is the analogue of Eq. (55) for the electroweak theory. For external Gold- 
stone fields extra contributions proportional to the Higgs vacuum expectation value 
appear, which are, however, irrelevant in the high-energy limit. Using (223), the LSC 
logarithms in (217) can be written as a single sum over external legs, 

gLSC ^h...i„ ^ ^ S^^^{k)Mo-"''-"\ (224) 



fc=i 



After evaluating the sum over A, Z, and W, in (217), the correction factors read 
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The first term represents the DL symmetric-electroweak part and is proportional to 
the electroweak Casimir operator C*^"" defined in (297). This is always diagonal in 
the SU(2) indices, except for external transverse neutral gauge bosons in the physical 
basis (301), where it gives rise to mixing between amplitudes involving photons and 
Z bosons. The second term originates from Z-boson loops, owing to the difference 
between Mw and Mz, and 

L-(s, A^ ml) := 2l(s) log (^^) + L(M^, A^) - L{ml A^) (226) 

contains all logarithms of pure electromagnetic origin. The LSC corrections for external 
longitudinal gauge bosons are directly obtained from (225) by using the quantum 
numbers of the corresponding Goldstone bosons. 

The contribution of the second term of (221) to (217) remains a sum over pairs of 
external legs, 

n 

Sssc^n.... ^ ^ ^ ^ ^ssc^j^^ l)Mi^-^-r-^" (227) 

fe=l Kk Va=A,Z,W^ 

with angular-dependent terms. The exchange of soft, neutral gauge bosons contributes 
with 

<SS^^'^) = 2 W)+l{M^^')] log^/4(A;)/4(/), 



5^^{k,l) = 2Z(.)log^i|,^(^)i|,(0 (228) 



and, except for in the neutral scalar sector H,x (see 6.1), the couplings are 
diagonal matrices. The exchange of charged gauge bosons yields 

= 2^(«)log^/|,(fc)i5,(0 (229) 

and owing to the non-diagonal matrices I^(k) [cf. (304), (309) and (313)], contributions 
of SU(2)-transformed Born matrix elements appear on the left-hand side of (227). In 
general, these transformed Born matrix elements are not related to the original Born 
matrix element and have to be evaluated explicitly. 

The SSC corrections for external longitudinal gauge bosons are obtained from (227) 
with the equivalence theorem (222) , i.e. the couplings and the Born matrix elements 
for Goldstone bosons^ have to be used on the right-hand side of (227). 

^Note that for Goldstone bosons x, the equivalence theorem as well as the couplings (310) and 
(308) contain the imaginary constant i. 
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The application of the above formulas is illustrated in Section 4 for the case of 
4-particle processes, where owing to ri2 = r^^, ri3 = and ru — ^23, (227) reduces 
to 

^SSC^ii«2i3«4 _ 
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and the logarithm with r^i = s vanishes. Note that this formula applies to 4 — * pro- 
cesses, where all particles or antiparticlcs and their momenta are incoming. Predictions 
for 2 2 processes are obtained by substituting outgoing particles (antiparticlcs) by 
the corresponding incoming antiparticlcs (particles). 

In addition to the angular terms in Eq. (230) at one looj) there are also SL terms 
from logarithms connected to parameter renormalization. These include the RG cor- 
rections to the mixing angles, electromagnetic charge, Yukawa couplings and scalar self 
couplings in the on-shell scheme. The corresponding terms are given in Ref. [29]. 

In Refs. [37, 38] we have shown that the splitting function approach employed for 
the electroweak radiative correction in the high energy regime described by Csymm is 
indeed justified at the one loop level. For this cross check it is crucial that the DL and 
the non- Yukawa SL terms factorize with respect to the same group factor. As can be 
seen by the explicit results in Eqs. (128) and (156), these terms are given by 
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for scalars. The correct reproduction of the explicit high energy approximation based 
on the physical fields in Ref. [97] constitutes a very strong check on the overall ap- 
proach. Since also the corrections from region II) in Fig. 1 are found to agree with 
Ref. [97] (including soft photon bremsstrahlung) , the matching condition at the weak 
scale is confirmed at the subleading level for both transverse as well as longitudinal 
degrees of freedom. The latter includes in particular the treatment of the would-be 
Goldstone bosons in the on-shell scheme to SL accuracy at high energies according to 
the effective Lagrangian Csymm in Eq. (11). 

This concludes our discussion of the one loop case. More definitions and represen- 
tations in the physical basis are presented in appendix 6.1. 
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Figure 20: The one loop clcctrowcak SM Fcynman diagrams leading to DL corrections 
in the 't Hooft-Feynman gauge for g — > /r/l- Only the vertex corrections from the 
neutral Z-boson (zigzag-lines) and the photon propagators contribute. At higher orders 
only corrections to these two diagrams need to be considered in the DL approximation. 
The photonic corrections are regulated by a fictitious mass terms A. In physical cross 
sections, the A-dependence is canceled by the eff'ect of the emission of soft and coUinear 
bremsstrahlung photons. 

3.9.2 Two loop results 

Also at the two loop level there have been checks performed with fixed order calcu- 
lations in terms of the physical SM fields. While the level of agreement demonstrated 
at the one loop level is impressive, it is very important for the desired accuracy of 
future linear colliders that there is agreement with real two loop DL calculations and 
the IREE method in order to trust predictions at the SL level. 

We follow Refs. [30, 36] to demonstrate the method used and to show how non- 
exponentiating terms cancel in the electroweak theory in the 't Hooft-Feynman gauge. 
Analogous results have been obtained in the Coulomb gauge in Ref. [31, 32] for external 
fermion lines and recently also for longitudinal and transverse gauge bosons. In this 
gauge, the DL corrections are related to two point functions and thus directly to 
external legs. The intermediate steps, however, are more cumbersome and therefore 
we only discuss the 't Hooft-Feynman gauge. In this gauge, a complete set of Feynman 
rules is listed in appendix 6.2. 

We begin with the case for right handed fcrmions coupled to a gluon, or in general a 
gauge singlet. At the one loop level, the electroweak corrections are depicted in Fig. 20. 
The fermion masses are neglected for simplicity. They can, however, be added without 
changing the nature of the higher order corrections. For right handed fermions we only 
need to consider the neutral electroweak gauge bosons, i.e. we are concerned with an 
?7(1)r X U{1)y gauge theory which is spontaneously broken to yield the Z-boson and 
photon fields. The DL-contribution of a particular Feynman diagram is thus given by 

Mk = MBorn (233) 

where the are given by integrals over the remaining Sudakov parameters at the 
n-loop level: 
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Figure 21: The two loop electroweak SM Feynman diagrams leading to DL corrections 
in the 't Hooft-Feynman gauge for g — > /r/l- The neutral Z-boson (zigzag-lines) and 
the photon propagators possess different on-shell regions due to the mass gap. 



The ©fe describe the regions of integration which lead to DL corrections. At one loop, 
the diagrams of Fig. 20 lead to 
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which is the well known result from QED plus the same term with a rescalcd coupling 
(see the Feynman rules in appendix 6.2) and infrared cutoff. The restriction to right 
handed fermions allows us to focus solely on the mass gap of the neutral electroweak 
gauge bosons. The only couples to left handed doublets. At the two loop level we 
have to consider more diagrams than in the QED case. The relevant Feynman graphs 
that give DL corrections in the 't Hooft-Feynman gauge are depicted in Fig. 21. Only 
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these corrections can yield four logarithms at the two loop level in the 't Hooft-Feynman 
gauge. Otherwise one cannot obtain the required pole terms (as is well known in QED 
[95]). They contain diagrams where the exchanged gauge bosons enter with differing 
on-shell regions, i.e. differing integration regions which give large DL corrections. It is 
instructive to revisit the case of pure QED corrections, since the topology of the graphs 
yielding DL contributions in the 't Hooft-Feynman gauge is unchanged. In QED at the 
two loop level, the scalar integrals corresponding to the first row of Fig. 21 are given 
by: 
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Thus, in QED we find the familiar result 
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whichs yields the second term of the exponentiated one loop result in Eq. (235) for 
Sw — > 0. In the electroweak theory, we also need to consider the remaining diagrams 
of Fig. 21. The only differences occur because of the rescaled coupling according to 
the Feynman rules in appendix 6.2 and the fact that the propagators have a different 
mass. Thus the second row of Fig. 21 leads to 
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where we indicate the gauge boson masses of the two propagators in the scalar func- 
tions. Analogously we find for the remaining two rows 

i>i' + S2 / / / OisuiVi - X )9{su2V2 - M )9(ui - U2) 

C4 ^0 Ui Jo Vi Jo U2 Jo V2 

^^Af du2 ^ifw2A_w2£!f£ 

Cl J Ml U2 JMl V2 2\ ^ X^ ^ A2 



80 



and 

rrMjM , r(M,M 
'-'l + '-'2 



1 s'^ 
2cl 



1, 2 



Thus, we find for tlie full two loop electroweak DL-corrections 
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which is precisely the second term of the exponentiated one loop result in the process 

of g — > fufi^- In Ref. [30] we showed that by using the appropriate quantum numbers 
listed in appendix 6.1, that the result of the IREE method of section 3.1 gives indeed 
the same result. 

For left handed fermions the calculation was performed in the 't Hooft-Feynman 
gauge in Ref. [36] . It is instructive to hst the intermediate steps of the calculation in 
order to identify terms not present in the case of unbroken gauge theories. We write 
the form of the correction as 



where the one loop result is given by 
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The two loop contribution can be written as a sum of the straight ladder diagrams, 
the crossed ladder and the three boson diagrams, which also enter due to the charged 
gauge boson exchange. Respectively we denote these terms as 
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The reason for this is that in QCD, the non-Abehan contribution of cl (~ Ca) is 
canceled by the tb terms. The exphcit results are 
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Thus we see that we have 
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which is precisely half the square of the first order correction. The novel feature here is 
that the way the cancellation happens is more complicated in the electroweak theory 
as already the straight ladder diagrams possess non-exponentiating terms. Only the 
full sum displays the familiar exponentiation. In the framework of the IREE method, 
however, these complications are naturally accounted for by carefully examining the 
kernel in both the high energy regime according to >Csymm in Eq- (H) and the region 
where only QED effects enter. At the scale M the appropriate matching conditions 
have to be employed. On the other hand, this complication is the reason for the wrong 
result in Ref. [25] which is based on QCD factorization properties. 

4 Applications 



In this section we apply the results presented in section 3 to specific processes at 
a future e+e" collider. As mentioned in the introduction, such a machine must be 
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a precision tool to disentangle and clarify the physics discovered at the LHC or the 

Tevatron. For our purposes we define precision as effects that change cross sections 
at the 1 % level. We would like to emphasize that at that level it is not sufficient 
to calculate only logarithmic corrections to all orders. At one loop, constant terms 

2 

such as ^ can lead to non-negligible effects and thus, a full one loop analysis is 
needed in these cases. In this work we want to emphasize the effect of the higher 
order corrections where constant terms are negligible. For all SL terms we can use the 
results summarized in Eq. (173) and the angular terms as described in section 3.6. In 
the following we will display only the "pure" electroweak corrections originating above 
the weak scale M relative to the Born cross sections. The importance of the angular 
dependent corrections is also discussed in Refs. [27, 110, 111] at higher orders and 
in [104] at the one loop level. While we focus below on the linear collider case, also 
for hadronic machines the electroweak corrections in the TeV range are important and 
must be included. In Ref. [119] the effect of one-loop logarithmic electroweak radiative 
corrections on WZ and Wj production processes at the LHC was studied for instance. 
Using the leading-pole approximation these corrections were implemented into Monte 
Carlo programs for pp luil'l', lui'-f. The authors find that electroweak corrections 
lower the predictions by 5-20 % in the physically interesting region of large transverse 
momentum and small rapidity separation of the gauge bosons. 

We begin in section 4.1 with transverse and longitudinal gauge boson production, 
followed by Higgs boson production in section 4.2 and heavy quark production in 4.3. 
For the latter two cases we also include supersymmetric Sudakov effects, in particular 
the dependence on tan/?. For all two to two processes we denote 

-^{l-cos9) , u = -^{l + cos9) (252) 

which is valid at high energies and where 9 is the scattering angle between the initial 
and final state particles. For definiteness, we use the following parameters 

mt = 174 GeV, M = 80.35 GeV, Z = 0.22356, — ^ ^ = (253) 

t , , w ' 47r 128 ^ ^ 

in all numerical results displayed. Furthermore we use g = e/s^ and g' = e/c^ below. 
The angular terms at one loop are obtained from Eq. (230) which was derived in Ref. 
[29] except for the case of Higgs-strahlung. All energy units in the figures are in GeV. 

4.1 Gauge boson production 

In this section we discuss the effect of the higher order pure electroweak corrections 
to transverse and longitudinal gauge bosons production at the cross section level. These 
processes, in particular the longitudinal ones, are important at high energies for the 
unitarity of the theory. If no Higgs boson should be found, the sector has to 
become strongly interacting in the TeV range in order to preserve unitarity. It would 
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Figure 22: The purely electroweak virtual corrections relative to the Born cross section 
in transverse and longitudinal production in e"'"e~ collisions at 90° scattering angle 
as a function of the cm. energy. The polarization is indicated in the figure for each 
symbol. Given are the one loop and the resummed corrections to SL accuracy in 
each case. It is clearly visible that the difference between the two approaches is non- 
negligible at TeV energies and necessitates the inclusion of the higher order terms. 
Pure QED corrections from below the weak scale are omitted. 
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therefore be important to know when and how the new dynamics would deviate from 
the perturbative SM predictions. 

Fig. 22 displays the energy dependence of the respective cross sections in 
production for the various polarizations. All universal terms are given by Eqs. (178), 
(179) and (177). In addition we have angular terms which read relative to the Born 
amplitude in the high energy effective regime described by Csymm in Eq. (11): 
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In terms of their numerical coefficients relative to the Born cross section we have 



da^'"^- = da^-^ -8.95 



t 



u 



log - + 1 - - log 



u 



-t 



""^Born 



S 

5.76 log - + 13.9 log — 



e2(M2), s ' 



u 



"•^Born 



-2.58 log 



t e\M' 



u 



87r2 M2 



dal'"'^'^ 

and RG corrections which at one loop are given by 
da 
da 
da 



87r2 



(257) 
(258) 
(259) 



eL,WT 
RG 






eL,WL 


— ""^Born ^ 


\ Al-82cl + 224 g\M') 


RG 


[ 124 47r2 


RG 


— «CrBoj.n 


f41 9'\M^) s 1 
^ 12 47r2 ^ M2 J 



(260) 

(261) 

(262) 
(263) 



The corrections are written relative to the Born cross sections and Fig. 22 shows the 
one loop and all orders results in each case to SL accuracy. The largest effect can be 
seen in the transverse sector for left handed electrons where the resummed terms are 
of order 30 % at 1 TeV, 55 % at 3 TeV and 65 % at 5 TeV. The respective cross section 
for right handed electrons is mass suppressed [97]. Also for longitudinal ( ~ 0^ 
via the equivalence theorem) we have significant effects. For left handed electrons the 
resummed terms are of order 20 % at 1 TeV, 38 % at 3 TeV and 46 % at 5 TeV. For 
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right handed electrons the resummed corrections are of order 11 % at 1 TeV, 22 % at 
3 TeV and 26 % at 5 TeV. 

In each case it can be seen that the one loop contributions alone are insufficient in 
the TeV regime and differ from the resummed results in the percentile range at 1 TeV, 
of order 10 % at 3 TeV and even more beyond. 

It is, however, not only important to take higher order DL corrections into account 
but also the higher order SL terms. This can be seen in Fig. 23 where the one loop 
terms agree (to SL accuracy) and at higher orders in one case only the DL and in the 
other also the SL terms are kept. Again the largest effect occurs in the transverse 
sector. The difference between the higher order SL and DL corrections is of order 4 % 
at 1 TeV, 10.5 % at 3 TeV and 14.5 % at 5 TeV. Also for longitudinal production 
the effects are significant. For left handed electrons the difference between the higher 
order SL and DL terms changes cross sections by 1.8 % at 1 TeV, 3 % at 3 TeV and 
3.5 % at 5 TeV. For right handed electrons we have differences of 0.5 % at 1 TeV, 2 % 
at 3 TeV and 3 % at 5 TeV. 

Thus, at these energies the SL terms can be as large as the leading terms and 
must not be omitted. For the longitudinal gauge bosons, there is a partial cancellation 
between the subleading Sudakov and Yukawa enhanced terms. Therefore, the overall 
effect is larger in the transverse sector but still large in the longitudinal one. 

Also the angular terms are significant as can be seen in Fig. 24. We treat the 
angular terms at higher orders as described in section 3.6 and the one loop terms from 
Eqs. (257), (258) and (259). The figure displays the effect for fixed cm. energy of 1 and 
3 TeV and we have written all angular terms in such a way that they are proportional 
to the Born cross section. This is always possible but involves factor of ^ etc. It can be 
seen that the angular corrections are large and vary for the resummed contributions by 
almost 20 % at 1 TeV and 17 % at 3 TeV for transverse production for scattering 
angles between 50° and 130". Also for longitudinal W"^ production the corrections are 
large. For left handed electrons the cross section changes by about 19 % at 1 TeV and 
20 % at 3 TeV in the same angular range. For the same range of scattering angles, the 
right handed electrons the resummed angular terms change cross sections by about 4.5 
% at 1 TeV and 5 % at 3 TeV. The one loop corrections are even larger and lead to 
significantly different results. Thus, the higher order terms are very important and it 
is mandatory to investigate if the two loop angular corrections in the full electroweak 
theory are indeed given by the product of the one loop terms and the Sudakov form 
factor. 

In e'^e^ collisions, the only non-mass suppressed longitudinal Z process is the 
Higgs-strahlung process. As such we discuss it in the next section together with other 
Higgs production processes. 
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Figure 23: The purely electroweak virtual corrections relative to the Born cross section 
in transverse and longitudinal production in e'^e" collisions at 90° scattering angle 
as a function of the cm. energy. The polarization is indicated in the figure for each 
symbol. Given are the full one loop and the resummed corrections to DL and SL accu- 
racy in each case. It is clearly visible that the difference between the two approaches, 
originating from SL terms at the two loop level, is non-negligible at TeV energies and 
necessitates the inclusion of the higher order SL terms. Pure QED corrections from 
below the weak scale are omitted. 
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Figure 24: The purely electroweak virtual corrections relative to the Born cross section 
in transverse and longitudinal production in e^e~ collisions at 1 and 3 TeV cm. 
energy as a function of the scattering angle. The polarization is indicated in the figure 
for each symbol. Given are the one loop and the resummed corrections to SL accuracy 
in each case. It is clearly visible that the difference between the two approaches is 
non-neghgible at TeV energies and necessitates the inclusion of the higher order terms. 
Pure QED corrections from below the weak scale are omitted. 
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4.2 Higgs production 



In this section we discuss light SM Higgs production processes relevant to the linear 
collider program. We begin with the Higgs-strahlung process e+e~ — Hx- The cross 
section for this process is smaller than the W^-fusion process at TeV energies, however, 
it would still be utilized to collect more Higgs events. The electroweak corrections, 
however, will reduce the cross section. 

All universal corrections are given by Eqs. (178) and (177). In addition we have 
angular terms which read relative to the Born amplitude in the high energy effective 
regime described by Csymm in Eq. (11): 
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Numerically we have for relative to the Born cross section: 
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and RG corrections which at one loop are given by 
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Fig. 25 depicts the changes of the cross section as a function of cm. energy for the 
two electron polarizations at 90° scattering angle. Only Z-cxchangc contributes. The 
resummed corrections reduce the cross section by about 11 % at 1 TeV, 23 % at 3 TeV 
and about 28 % at 5 TeV for left handed electrons. For right handed electrons, the 
effect is actually an increase in the cross section at energies up to about 2 TeV due to 
the large positive angular terms. The reduction at 3 TeV is about 5 % and 14 % at 5 
TeV for the resummed cross sections. Again wc can sec that the one loop predictions 
differ in the percentile regime at 1 TeV and by about 7 % at 3 TeV. At higher energies 
the difference grows even more rapidly. Thus, for the linear collider in the TeV range, 
higher order contributions are necessary. 

In Fig. 26 the angular dependence for the Higgs-strahlung process is displayed. 
Only the cross section involving left handed electrons possess angular dependent terms. 
The corrections are symmetric with respect to the central scattering angle (i.e. sym- 
metric in u ^ t). The resummed cross sections are consistently lower due to the large 
positive angular one loop result in Eq. (265). At 1 TeV, the corrections change by 
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Figure 25: The purely electroweak virtual corrections relative to the Born cross section 
in Higgs-strahlung in e+e" collisions at 90° scattering angle as a function of the cm. 
energy. The polarization is indicated in the figure for each symbol. Given are the one 
loop and the resummed corrections to SL accuracy in each case. The difference between 
the two approaches is non-negligible at TeV energies and necessitates the inclusion of 
the higher order terms. Pure QED corrections from below the weak scale are omitted. 
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Figure 26: The purely electroweak virtual corrections relative to the Born cross section 
in Higgs-strahlung in e^e~ collisions at 1 and 3 TeV cm. energy as a function of the 
scattering angle. Only cross sections involving left handed electrons receive SL angular 
corrections. Given are the one loop and the resummed corrections to SL accuracy in 
each case. The difference between the two approaches is non-negligible at TeV energies 
and necessitates the inclusion of the higher order terms. Pure QED corrections from 
below the weak scale are omitted. 
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Figure 27: The purely electroweak virtual corrections relative to the Born cross section 
in Higgs radiation off (fusion) in e+e" collisions as a function of the cm. energy. 
The polarization is indicated in the figure for each symbol. Given are the one loop and 
the resummed corrections to DL accuracy in each case. It is clearly visible that the 
difference between the two approaches is non-negligible at TeV energies and necessi- 
tates the inclusion of the higher order terms. Pure QED corrections from below the 
weak scale are omitted. The SL terms are omitted for simplicity, however, are needed 
for a full treatment. 
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about 14 % and at 3 TeV by about 18 % over the displayed angular range. Again we 
sec how important those terms are for future collider experiments in the TeV range. 

We now want to briefly discuss two other important Higgs production processes 
at e'^e~ colliders, namely the W fusion process and Higgs-strahlung off top quarks. 
Both of these processes involve three final state particles and given the multiplicity 
of final state scattering states, we restrict ourselves here to a discussion only of the 
energy-dependence at DL accuracy. All universal DL corrections are given by the DL 
terms in Eqs. (128) and (156). 

We begin with the l^-fusion process e~^e~ — > i^e^eH in Fig. 27. Only left handed 
electrons contribute. This process yields the largest cross section at high energies for 
Higgs production. The resummed DL corrections reduce the cross section by about 25 
% at 1 TeV, 46 % at 3 TeV and 53 % at 5 TeV. The one loop DL predictions are off 
by about 5 % at 1 TeV, 15 % at 3 TeV and 20 % at 5 TeV relative to the Born cross 
section! 

While the SL terms can change this picture significantly, the typical size of the 
corrections should be correctly depicted in Fig. 27. 

In Fig. 28 we display the electroweak DL corrections to Higgs productions off top 
quarks for the various polarizations. This process is crucial in order to determine the 
top-Higgs Yukawa coupling [120]. The DL terms reduce the cross section for cl^l by 
about 24 % at 1 TeV, 43 % at 3 TeV and 52 % at 5 TeV. For CRtL by about 18 % at 1 
TeV, 35 % at 3 TeV and 42 % at 5 TeV. For CLtR by about 17 % at 1 TeV, 33 % at 3 
TeV and 40 % at 5 TeV and finally for CRtR by about 11 % at 1 TeV, 21 % at 3 TeV 
and 27 % at 5 TeV. Again the one loop DL corrections differ significantly, especially 
for left handed polarizations. 

In all Higgs production processes discussed in this section, the electroweak radiative 
corrections are important and can reduce the cross sections considerably. Even for 
the top- Yukawa measurement at 800 GeV at TESLA or the NLC, corrections are of 
O (20%) and the difference between one loop and resummed and be a few percent. 
Therefore a full higher order SL analysis is warranted for this process in addition to 
the QCD corrections. 

4.2.1 Charged MSSM Higgs production 

In this section we discuss the effect of MSSM Sudakov effects in charged Higgs 
production. The relevant radiative corrections are given in Eq. (199) for the universal 
and angular dependent Sudakov terms above the susy scale set by = tuh under 
the assumptions stated in section 3.7. The angular dependent corrections can also 
be obtained form Eqs. (258) and (259) via the replacement M — > tjih, which for our 
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Figure 28: The purely electroweak virtual corrections relative to the Born cross section 
in Higgs radiation off final state top-quarks in e+e" collisions as a function of the cm. 
energy. The polarization is indicated in the figure for each symbol. Given are the one 
loop and the resummed corrections to DL accuracy in each case. It is clearly visible 
that the difference between the two approaches is non-negligible at TeV energies and 
necessitates the inclusion of the higher order terms. Pure QED corrections from below 
the weak scale are omitted. The SL terms are omitted for simplicity, however, are 
needed for a full treatment. 
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purposes, however, is of SSL accuracy in our "light susy" mass assumption. In addition 
we have RG contributions at one loop which read 

' " I n{a{mj,) + a'{mj,)) j^''^^ ^^^^^ 

daT'^- = dal^Z' [-2^^i3'}> log ^ (270) 

where /3o and {3^ are given in Eq. (198). In Fig. 29 the energy dependence for the 
process e^e~ — > H^H~ is depicted for two typical values of tan /5 = 10 and tan {3 = 40. 
The scattering angle is held fixed at ^ = 90" since the angular dependence is analogous 
to the case of longitudinal W production in Fig. 24. The helicity of the electrons is 
indicated and only the resummed all orders SL results are presented. The value of the 
heavy charged Higgs particles is rriH = 300 GeV. For a different heavy Higgs mass, the 
displayed results shift accordingly. 

Fig. 29 shows that the effect of the MSSM Sudakov corrections is somewhat 
reduced compared to the SM case. This is almost entirely due to the different mass 
scale used. The left handed cross sections are more suppressed than the right handed 
ones due to the larger DL group factors. In addition, larger values for tan/3 lead to 
an enhanced suppression since the overall sign of the Yukawa terms is negative. The 
dependence on tan /3 is explicitly shown in Fig. 30 for cm. energies of a/s — 1 TeV 
and ^/s — 3 TeV. At 1 TeV the right and left handed cross sections display the same 
overall dependence on tan/3, differing by about 7 % in the range displayed. At 3 TeV, 
the left handed cross section is more suppressed but the overall tan/5 dependence is 
similar. The cross sections can differ by about 13 % in the range between tan/3 = 6 
to tan (3 — 40. This strong dependence, especially for larger values on tan /3 originates 
from both, the size of the Yukawa terms as well as from the number of colors in the 
quark loops. In Ref. [114] this was utilized to suggest a measurement of tan/3 at CLIC 
with a relative precision better than 25 % (a few percent for large values tan/3 > 25). 
More importantly, this determination of tan /3 to SL accuracy does not depend on the 
soft breaking terms (which are constants) and is scheme and gauge invariant. 

While we focussed here only on the case of heavy charged Higgs production, also the 
neutral Higgs, fermion and sfermion processes have important information on MSSM 
parameters contained in the coefficients of large Sudakov logarithms and should be 
fully exploited at such a collider. 

4.3 Heavy quEirk production 

In this section we discuss heavy quark production at the linear collider. These 
processes can be used to measure as above the production threshold and should be 
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Figure 29: The purely electroweak MSSM virtual corrections to charged Higgs produc- 
tion above the susy scale set by inn = 300 GeV to SL accuracy. The dependence in 
shown for two characteristic choices of tan j3 for both right and left handed electron 
polarizations. 
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Figure 30: The purely electroweak MSSM virtual corrections for charged Higgs pro- 
duction above the susy scale set by itlh = 300 GeV to SL accuracy as a function of 
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fully understood. A general result for the angular terms relative to the Born amplitude 
is given by 
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log ^ \ [tan^ e^Y^-Y,, + 4r^_r|J log i 



where the last line only contributes for left handed fermions and the d, u symbols 
denote the corresponding isospin quantum number of /. 

We begin with top production. All universal terms are given by Eqs. (178) and 
(177). The angular terms and RG corrections which at one loop are given by 
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Fig. 31 displays the energy dependence of the corrections for central scattering angles 
for the various polarizations. The difference between the resummed and the one loop 
contributions is in the several percent range only above 1 TeV. The largest corrections 
are again obtained in the case where both fermions are left handed. The overall 
corrections are large and non- negligible. 

Also the angular terms, depicted in Fig. 32, arc significant, especially for eiti 
where they change cross sections by about 20 % for 1 and 3 TeV over the angular 
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Figure 31: The purely electroweak virtual corrections relative to the Born cross section 
in top quark production in e+e" collisions at 90° scattering angle as a function of the 
cm. energy. The polarization is indicated in the figure for each symbol. Given are the 
one loop and the resummed corrections to SL accuracy in each case. The difference 
between the two approaches is non-negligible at TeV energies and necessitates the 
inclusion of the higher order terms. Pure QED corrections from below the weak scale 
are omitted. 
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Figure 32: The purely electroweak virtual corrections relative to the Born cross section 
in top quark production in e'^e" collisions at 1 and 3 TeV cm. energy as a function of 
the scattering angle. The polarization is indicated in the figure for each symbol. Given 
are the one loop and the resummed corrections to SL accuracy in each case. The dif- 
ference between the two approaches is non-negligible at TeV energies and necessitates 
the inclusion of the higher order terms. Pure QED corrections from below the weak 
scale are omitted. 
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Figure 33: The purely electroweak virtual corrections relative to the Born cross section 
in bottom quark production in e+e" collisions at 90" scattering angle as a function of 
the cm. energy. The electron polarization is indicated in the figure for each symbol. 
Given are the one loop and the resummed corrections to SL accuracy in each case. 
The difference between the two approaches is non-negligible at TeV energies and ne- 
cessitates the inclusion of the higher order terms. Pure QED corrections from below 
the weak scale are omitted. 
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range between 50" and 130°. This is another indication that the higher order angular 
terms need to be understood in the full electroweak theory. 

For bottom quark production does not contain SL- Yukawa terms. All universal 
corrections are given by Eqs. (178) and (177). In addition we have angular terms and 
RG corrections which at one loop are given by 



da, 



da^"''"' 



da' 



.eL,bL 
Born 



4.04 log - + 19.8 log — 

u s 



e\M' 



87r2 M2 



= d<7^-^^ 0.86 log - 



t e\M^) 



^o^a^omN -0.86 log 



u 87r2 M2 ^ 

t e\M^) 



""^Born 



U 



87r2 M2 



-1.72 log- 



t e\M^) 



u 87r2 M2 



(280) 
(281) 
(282) 
(283) 



and 



da 



da 



da 



eL,bL 
RG — 


""^Born ^ 


e2(M2) 


RG — 


""^Born 




eL,bR 
RG — 


""^Born 




eR,bR _ 
RG — 


""^Born 






(284) 
(285) 
(286) 
(287) 



Fig. 33 demonstrates that in the energy range displayed, the electroweak corrections 
are actually positive for enhR. This is mainly due to the RG corrections for the right 
handed coupling {g') which is Abelian and therefore increases with energy. Secondly 
it is due to the fact that the Yukawa terms are absent and the DL terms are partially 
offset by the SL terms. 

The remaining corrections are similar to the top production discussion. 

The angular terms, depicted in Fig. 34, show a similar behavior to the case of top 
quark production in Fig. 32 in that the largest contribution is again for the purely left 
handed case. The overall size of the corrections differs and the euhu contributions stay 
positive for most of the displayed angular range. 

5 Outlook 



Electroweak radiative corrections at high energies have received much attention 
recently due to their importance at experiments in the TeV regime. It is not only the 
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(da^™-daB)/claB 
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A f 



; e^^b^ exp () 
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: Cjb^ lloop 
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□ : Bj^g lloop 

A : Bjbg exp () 

: ejfg lloop 
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Figure 34: The purely electroweak virtual corrections relative to the Born cross section 
in bottom quark production in e^e~ collisions at 1 and 3 TeV cm. energy as a function 
of the scattering angle. The electron polarization is indicated in the figure for each 
symbol. Given are the one loop and the resummed corrections to SL accuracy in each 
case. The difference between the two approaches is non- negligible at TeV energies and 
necessitates the inclusion of the higher order terms. Pure QED corrections from below 
the weak scale are omitted. 
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phenomenological importance, however, that has led to a surge in interest into the high 
energy behavior of the SM but also the fact that conceptually new effects enter due to 
EWSB. The main differences to the case of unbroken gauge theories originate from the 
longitudinal sector and the fact that three of the gauge bosons acquire masses, while 
the photon stays massless. In addition, the asymptotic states carry a non-Abelian 
group charge, the weak isospin and are superpositions of the fields of the unbroken 
phase. 

As a consequence, fully and semi-inclusive cross sections show double and single 
logarithmic dependencies on the ratio of the energy and the gauge boson mass, and 
longitudinal degrees of freedom are not mass suppressed. Thus, observables in the SM 
depend on the infrared cutoff, the weak scale, which in this case, however, is a physical 
parameter. At one loop, general methods exist which allow to calculate corrections 
relevant to precision measurements at future colliders. Also at higher orders, a general 
approach to SL accuracy is available, based on the high energy approximation of the 
symmetric part of the SM Lagrangian and is phenomenologically necessary. 

Focusing on techniques to calculate the higher order corrections, we have summa- 
rized the present status of virtual electroweak radiative corrections to SM and MSSM 
high energy processes. In the framework of the IREE method, the high energy effective 
theory is based on the high energy limit of the SM Lagrangian in the symmetric limit 
where all terms connected to the v.e.v. can be neglected to SL accuracy. The QED 
corrections below the weak scale are incorporated with the appropriate matching con- 
ditions. This approach is so far the only one able to allow for a two loop calculation 
of the DL and SL terms which are relevant for TeV experiments at future colliders. 
We have shown that the one loop terms are insufficient when the cm. energy is larger 
than 1 TeV and that both, DL and SL terms at the two loop level are necessary at 
the several percent level. Also angular dependent corrections cannot be neglected and 
it should be investigated if their calculation at the SSL level is needed for some ob- 
servables. It should also be investigated if large Yukawa constants can be treated in 
a systematic manner to SSL accuracy. In this context it is also important to consider 
the emission of real gauge bosons above the weak scale even at the SL level. 

In summary, there exists a way to calculate all higher order virtual SL electroweak 
radiative corrections to high energy processes. The approach is in agreement with all 
available one loop calculations in terms of the physical SM or MSSM fields and at the 
two loop level to DL accuracy. These terms are crucial for the experiments at future 
colliders in the TeV regime since the effects of new physics expected in this range can 
be rather small. 
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6 Appendix 



In this appendix we list the relevant quantum numbers of the physical SM fields 
and the corresponding Feynman rules of the full SM in the 't Hooft-Feynman gauge. 

6.1 Operators and quantum numbers 

In this section we present the generators of the physical gauge group and various 
group-theoretical matrices used in section 3.9. We follow the discussion of Ref. [29]. 
The notation for the components of such matrices is 

M^,^M (288) 

where the argument ip represents a multiplct and fixes the representation for the matrix 
M, whereas 99, are the components of the niultiplet. Explicit representations for left- 
and right-handed fermions {^p — , f^, , f^)-, for gauge bosons (</? = V) and for 
the scalar doublet [(p = $) arc given below. Where the representation is implicit, the 
argument (f is omitted. For the eigenvalues of diagonal matrices one has 

M^,^, =5^,^,M^^ (289) 

Symmetric and physical gauge fields and gauge couplings 

For physical gauge bosons one needs to take special care of the effect of Weinberg 
rotation (mixing). The symmetric basis Va = B,W^,W^,W^, is formed by the U(l) 
and SU(2) gauge bosons, which transform as a singlet and a triplet, respectively, and 
quantities in this basis are denoted by a tilde. The physical basis is given by the charge 
and mass eigenstates K = A, Z, W^, W". The physical charged gauge bosons, 

= , (290) 

are pure SU(2) states, whereas in the neutral sector the SU(2) and U(l) components 
mix, and the physical fields N = A, Z are related to the symmetric fields N — B, 
by the Weinberg rotation. 



N = Uj^MN, U{e^) ={ ^ " (291) 




with Cy, = cos6'„ and = sin6'w In the on shell renormalization scheme the Weinberg 
angle is fixed by 
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The gauge couplings are given by the generators of global gauge transformations. 
In the symmetric basis, they read 



a = 1,2,3 



(293) 



where Y is the weak hypercharge and are the components of the weak isospin. In 
the physical basis one has 



V2 



(294) 



with Q = T^ + Y/2. 

Casimir operators 

The SU(2) Casimir operator is defined by 



c = ^(r«)2 

a=l 



(295) 



Loops involving charged gauge bosons are often associated with the product of the 
non-Abelian charges 



{I 



W\2 



a=± 



C -{T 



3\2 



(296) 



and if one includes the contributions of neutral gauge bosons, one obtains the effective 
electroweak Casimir operator 



^ew ._ jVajVa 



Va=A,Z,W± 



(297) 



For irreducible representations (fermions and scalars) with isospin T^, the SU(2) 
Casimir operator is proportional to the identity and reads 



C^iip.,{<^) — S^^^p^,C^, — T^[T^p + 1] 



(298) 



Physical gauge bosons have a reducible representation as already discussed in section 
3. In the symmetric basis C{V) is a diagonal 4x4 matrix 



(299) 



with U(l) and SU(2) eigenvalues 



Cb — 0, 



(300) 
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The transformation of a matrix like (299) to the physical basis, yields a 4 x 4 matrix 
with diagonal 2x2 block structure, i.e. without mixing between the charged sector 
(W^^^) and the neutral sector {N = A,Z). In the neutral sector C{V) becomes non- 
diagonal owing to mixing of the U(l) and SU(2) eigenvalues. 



Cnn' — 



U{0^)CU-\9^) =2{ ';^| (301) 



whereas in the charged sector it remains diagonal, 

C^^a^ra' — ^.Scra' (302) 

Explicit values for Y, Q, T\ C, {I^f, {I^f, C^^, and /± 

Here a list of the eigenvalues (or components) of the operators F, Q, T^, C, (/^)^, 
(7^)^, (/^)^, C*^™, and is given following Ref. [29], that have to be inserted in 
the general results. For incoming particles or outgoing antiparticles the values for the 
particles have to be used, for incoming antiparticles or outgoing particles the values of 
the antiparticles. 



Fermions 

The fermionic doublets = {fl-fl)"^ transform according to the fundamental 
or trivial representations, depending on the chirality k = L, R. Except for I^, the 
above operators are diagonal. For lepton and quark doublets, L'^ = (i/**, l'^)'^ and 
QK ^ ^^K^ d'^Y , their eigenvalues are 



r/2 g r3 c {i^f {i^f {i^f c 










±1 


3 





1 


1 


1+24 


4 


4s^ 


2^ 






T| 


Tl 




3 
4 


1 


(c"^ _s2 -12 


1 
24 


1+24 






Tl 








1 


s2 
£SL 





1 






±1 


±1 


3 


4 


(3c^— sj,)^ 


1 




1 


9 






364 6-2, 


d^,d^ 


H 






3 


1 




1 




4 


9 




24 






±i 


±i 








4 
9 


4 4 
9 c2 





4 
9c^ 














1 

9 


1 si 





1 

9c2 



(303) 



For left-handed fermions, /^(/^) have the non- vanishing components 

'If-jn = = ^g; (304) 

whereas for right-handed fermions /^(/^) = 0. 
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Scalar fields 



The symmetric scalar doublet, $ = (0"^,0o)^, = (0~;0o)^' transforms accord- 
ing to the fundamental representation, and its quantum numbers correspond to those 
of left-handed leptons (303) with 



(f)'^ l^, 00 2^^, 



0- ^ l^, 



00 



(305) 



After symmetry breaking the neutral scalar fields are parametrized by the mass eigen- 
states 

cf>o=^{v + H + ix) (306) 

With respect to this basis S = {H,x) the operators Q,C,{I^)'^, and C^^ remain 
unchanged, while and Y become non-diagonal in the neutral components 



T, 



SS' 



I) = 

^ / SS' 



1 / -i 

2 i 



and 

The W='= couphngs read 
with 



I'' 
^S4>- 



Saa'Is 



2s. 



X 



2s. 



(307) 

(308) 

(309) 
(310) 



Gauge fields 

For transversely polarized external gauge bosons one has to use the adjoint repre- 
sentation. In the symmetric basis the diagonal operators have eigenvalues 





Y/2 


Q 




c 















±1 


±1 


2 


1 


i-w 


1 

s2 
*w 


2 
*w 













2 








2 


2 


B 



























(311) 



In the neutral sector, owing to the Weinberg rotation, the non-trivial operators C^^, C 
and [I^Y become non-diagonal in the physical basis N — A,Z, with components 



/^ew ^rl^^2 

^NN' — -T^NN' — \J- )nN' 



(312) 
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whereas the trivial operators Y/2 — Q — — {I'^Y — {I^Y = remain unchanged. 
In the physical basis the non- vanishing components of the couplings are 

^'nW-'^' ~ ~^W<^' N ~ ^(Tcr'^N (313) 

with 

II = Iz = ^- (314) 

6.2 Elect roweak Feynman rules 

In this appendix we list the Feynman rules of the SM in the 't Hooft- Feynman gauge 
including the counterterms in a way appropriate for the concept of generic diagrams 
[41]. I.e. we write down generic Feynman rules obtained from the classic Lagrangian 
in Eq. (5) and give the possible actual insertions. We omit any field rcnormalization 
constants for the unphysical fields. For brevity we introduce the shorthand notation 

c = Cw, s = Sw (315) 

In the vertices all momenta are considered as incoming. 
Propagators: 

for gauge bosons F = 7, Z, W in the 't Hooft-Feynman gauge (^j = 1) 











A;2 




for Faddeev-Popov ghosts G 














k 


-G 




i 

-Ml 


for scalar fields S = H, x, <P 












S — 
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— S 


k^ 


i 

-Ml 


and for fermion fields F — fi 














P 




i(^ 


+ mp) 










p2 





In the 't Hooft-Feynman gauge we have the following relations: 

Mui = 0, M„z = M^^Mz, M„± = = Mw (316) 
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Tadpole: 



VV counterterm: 



iSt. 




with the actual values of Vi, V2 and Ci, C2 

without renormalization of the gauge-fixing term 



w+w- 


: Ci 


-^3 = 




zz 


: Ci 


= C3 = 




AZ 


: Ci 


= C3 = 




AA 


: Ci 


= ^3 = 





C2 = M2,5ZW + 

C2 = M^SZ^^ + SMl 



C,. = Q 



with renormalization of the ('t Hooft-Feynman) gauge- fixing term 

C3 = 0, Ci, C2 as above 



VS counterterm: 




with the actual values of V, S and C 

without renormalization of the gauge-fixing term 

W±0T : C = ±Mw|(5ZW + ^), 
Zx : C = iMzi(<5Z^^ + ^), 
Ax : C = iMz^SZ^^ 

with renormalization of the ('t Hooft-Feynman) gauge- fixing term 

C = 
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SS-counterterm: 




CiA;2 - C2 



with the actual values of Si, S2 and Ci, C2 

without renormalization of the gauge-fixing term 

HH : Ci = SZh, C2 = Ml5Z^ + 5Ml, 

: Ci = 0, C2 = -^^ 
with renormalization of the ('t Hooft-Feynman) gauge- fixing term 
HH : Ci = 6Zh, C2 = M^6Z^ + SM^, 
XX ■■ Ci = 0, C2^-t^ + 5Ml 



FF-counterterm: 




with the actual values of Fi, F2 and Cl, Cr, Cg , Cg 



fjf^ 
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VWV coupling: 



> 


< 







with the actual values of VI, V2, V3, V4 and C 
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(321) 



VW couphng: 



= -ieC 



with the actual values of 14, V2, and C 
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(322) 
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SSSS coupling: 




with the actual values of Si, 82, S^, and C 
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(323) 



SSS coupling: 




ieC 



with the actual values of Si, S2, S3 and C 
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WSS coupling: 
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with the actual values of Vi, V2, Si, S2 and C 
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VSS coupling: 





Si, ki 
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with the actual values of V, Si , S2 and C 
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SW coupling: 
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VFF coupling: 
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(328) 
(329) 



The vector and axial vector couplings of the Z-boson are given by 

- U9J + 9t) - «/ = 1(^7 - 91) = % 



(330) 
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SFF coupling: 




ie(C-a;_ + C+uj+) 



with the actual values of S, Fi , F2 and C"*", C 
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VGG coupling: 




with the actual values olV, Gi , G2 and C 
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